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Abstract. We study an integrable vertex model with a periodic boundary condition associa- 
ted with Uqi^An^) at the crystallizing point q = 0. It is an (n + l)-state cellular automaton 
describing the factorized scattering of solitons. The dynamics originates in the commuting 
family of fusion transfer matrices and generalizes the ultradiscrete Toda/KP flow corre- 
sponding to the periodic box-ball system. Combining Bethe ansatz and crystal theory in 
quantum group, we develop an inverse scattering/spectral formalism and solve the initial 
value problem based on several conjectures. The action-angle variables are constructed rep- 
resenting the amplitudes and phases of solitons. By the direct and inverse scattering maps, 
separation of variables into solitons is achieved and nonlinear dynamics is transformed into 
a straight motion on a tropical analogue of the Jacobi variety. We decompose the level set 
into connected components under the commuting family of time evolutions and identify each 
of them with the set of integer points on a torus. The weight multiplicity formula derived 
from the q = Bethe equation acquires an elegant interpretation as the volume of the phase 
space expressed by the size and multiplicity of these tori. The dynamical period is deter- 
mined as an explicit arithmetical function of the n-tuple of Young diagrams specifying the 
level set. The inverse map, i.e., tropical Jacobi inversion is expressed in terms of a tropical 
Riemann theta function associated with the Bethe ansatz data. As an application, time 
average of some local variable is calculated. 

Key words: soliton cellular automaton; crystal basis; combinatorial Bethe ansatz; inverse 
scattering/spectral method; tropical Riemann theta function 
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1 Introduction 

1.1 Background 

In [25], a class of soliton cellular automata (SCA) on a periodic one dimensional lattice was 
introduced associated with the non-exceptional quantum affine algebras Ug{Qn) at q = 0. In 

this paper we focus on the type 0„ = An '^ . The case n = 1 was introduced earlier as the 
periodic box-ball system [41]. It is a periodic version of Takahashi-Satsuma's soliton cellular 
automaton [38] originally defined on the infinite lattice. In the pioneering work [iO], a closed 
formula for the dynamical period of the periodic box-ball system was found for the states with 
a trivial internal symmetry. 

After some while, the first complete solution of the initial value problem and its explicit 
formula by a tropical analogue of Riemann theta function were obtained in [241 [20l [2T] . This 
was done by developing the inverse scattering/spectral methocll] [6] in a tropical (ultradiscrete) 
setting primarily based on the quantum integrability of the system. It provides a linearization 
scheme of the nonlinear dynamics into straight motions. The outcome is practical as well 
as conceptual. For instance practically, a closed formula of the exact dynamical period for 



^Although it is more customary to use 'spectral' in periodic systems, we shall simply call it inverse scattering 
hereafter for short. 
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arbitrary states [24^ Theorem 4.9] follows straightforwardly without a cumbersome combinatorial 
argument, see (|3.4U|) . Perhaps what is more important conceptually is, it manifested a decent 
tropical analogue of the theory of quasi-periodic soliton solutions [3l 0] and its quantum aspects 
in the realm of SCA. In fact, several essential tools in classical integrable systems have found 
their quantum counterparts at a combinatorial level. Here is a table of rough correspondence. 



classical 


quantum 


solitons 


Bethe strings 


discrete KP/Toda flow 


fusion transfer matrix 


action-angle variable 


extended rigged configuration mod {q = Bethe eq.) 


Abel-Jacobi map, Jacobi inversion 


modified Kerov-Kirillov-Reshetikhin bijection 


Riemann theta function 


charge of extended rigged configuration 



Here the Kerov-Kirillov-Reshetikhin (KKR) bijection |15l [T6] is the celebrated algorithm on 
rigged configurations in the combinatorial Bethe ansatz explained in Appendix O As for the 
last line, see (jl.ip . 

In a more recent development [36], a complete decomposition of the phase space of the periodic 
box-ball system has been accomplished into connected components under the time evolutions, 
and each of them is identified with the set of integer points on a certain torus TP jF^L^ explicitly. 



1.2 A quick exposition 

The aim of this paper is to extend all these results |24l [20l [211 [36] to a higher rank case based 
on several conjectures. The results will be demonstrated with a number of examples. We call 
the system the •periodic soliton cellular automaton (SCA). It is a dynamical system on 
the L numbers {1,2,... , n + 1}^ which we call paths. There is the n-tuple of commuting series 
of time evolutions . . . , T^"^ with / > 1. The system is periodic in that the cyclic shift ()2.3p 

is contained in the member as T^^\ The casual exposition in the following example may be 
helpful to grasp the idea quickly. 

Example 1.1. Take the path p = 321113211222111223331111 of length L = 24 with n = 2, 
which will further be treated in Examples 12. 3[ \'6.'2\ 13.4113.51 and 14. 3[ 







{tP)\p) 


(Tf)*(p) 


t = 


: 


321113211222111223331111 


321113211222111223331111 


t = 


1 : 


113211132111222112213331 


132111321333111221112221 


t = 


2 : 


331132111321111221122213 


112211132111222331113321 


t = 


3 : 


213311332113211112211122 


213311113211322112211132 


t = 


4 : 


122233111332132111122111 


221122111321133213311112 


t = 


5 : 


111122332111321332111221 


331132111122111321122213 


t = 


6 : 


211111221332113211332112 


112213221133111132132211 


t = 


7 : 


122211112211332132111331 


113311222111221113213321 


t = 


8 : 


311122211122111321332113 


211122333111321111221132 


t = 


9 : 


133111122211221113211332 


221132111222132111331113 



By regarding 1 as an empty background, these patterns show the repeated collisions of solitons 
like 22, 32, 33, 222, 322, 332 and 333. At each time step there are always 5 solitons with 
amplitudes 3, 3, 2, 2, 2. The dynamics is highly nonlinear. Ultimately in our approach, it will 
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be transformed into the straight motion (Hint: right hand side of (j3.15p and Remark 13. ip : 



(r/^))*(T(?))^(p) 











1 + min(Z,3)i 


















+ min(Z,3)t 




1 + min 








7 + min(/,2)t 








4 + min(/,2)t 








2- 


f min(/,2)t 



11 + min(m, 4)s 



Here $ is the modified Kerov-Kirillov-Reshetikhin (KKR) bijection mentioned in the table. 
One notices that Tg^^ = and T^^^ = in this example. (Of course min(m, 1) = 1.) 
Note also that the first Young diagram (33222) gives the list of amplitudes of solitons. The two 
Young diagrams are conserved quantities (action variable). A part of a Young diagram having 
equal width is called a block. There are 4 blocks containing 2, 3, 1, 1 rows in this example. The 
7 dimensional vector consisting of the listed numbers is the angle variable flowing linearly. It 
should be understood as an element of (Hint: ()4.17p ) 

z7-Bz7(e2 X 63 X ©1 X ©1), 

where &m is the symmetric group of degree m. We have exhibited the trivial ©I's as another 
hint. The matrix B is given by (j4.1ip . The description is simplified further if one realizes that 
the relative (i.e., differences of) components within each block remain unchanged throughout. 
Thus picking the bottom component only from each block, we obtain 



( + min(/,3)t \ 
2 + min(/,2)i 
11 + min(m, 4)s 
\ 1 + min(m, l)s / 



mod Fl} 



F 



/16 12 -3 

8 19 -2 

-6 -6 10 

V-2 -3 2 



-1\ 

-1 

2 

3/ 



The 4 dimensional vector here is the reduced angle variable (Section l3.3p . After all, the dynamics 
is just a straight motion in Z^/i^'Z^. The abelian group of the time evolutions acts transitively 
on this torus. Thus the size of the connected component of p is det F = 4656. The dynamical 
period A/'/''^ of tI^^\ i.e., the smallest positive integer satisfying (T^*-''^)-^; (p) = p is just the 
smallest positive integer satisfying 



(1) 



/ min(/, 3)\ 
min(/, 2 





e F2 








min(?7T,, 4) 
\min(m, 1) / 



where the coefficient vectors are velocities. Solving the elementary exercise, one finds 





AA,« 


AAi^) 


AAf) 


AAi^) 


MP 




24 


12 


194 


1164 


776 


582 


2328 



The level set characterized by the above pair of Young diagrams ((33222), (41)) consists of 139680 
paths. It is decomposed into 36 connected components as (see Example 13. 6p 



24(Z7ra^) U 12(Z^/F'Z^), 139680 = 24(detF) + 12(detF'), 

where F' = F ■ diag(^, 1, 1, 1). Here the two kinds of tori arise reflecting that the two kinds of 
internal symmetries are allowed for the pair ((33222), (41)) under consideration. 
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1.3 Related works 

A few remarks are in order on related works. 

1. The solution of initial value problem of the n = 1 periodic box-ball system [241 [20l [2T] 
has been reproduced partially by the procedure called 10-elimination [29]. By now the precise 
relation between the two approaches has been shown [17j . It is not known whether the 10- 
elimination admits a decent generalization. The KKR bijection on the other hand is a canonical 
algorithm allowing generalizations not only to type An^ [18] but also beyond [32]. Conjecturally 
the approach based on the KKR type bijection will work universally for the periodic SCA |25j as- 
sociated with Kirillov-Reshetikhin crystals as exemplified for A^^^ higher spin case [22] and A^n^ 
[25 \ 1261 this paper]. In fact, the essential results like torus, dynamical period and phase space 
volume formula are all presented in a universal form by the data from Bethe ansatz and the 
root systenl^. From a practical point of view, it should also be recognized that the KKR map 
(f>^^ is a delightfully elementary algorithm described in less than half a page in our setting in 
Appendix [O It is a useful exercise to program it to follow examples in this paper. 

2. When there is no duplication of the amplitudes of solitons, the tropical analogue of 
the Jacobian i7(/u) obtained in [24] has an interpretation from the tropical geometry point of 
view [9] |30]. In this paper we are naturally led to a higher rank version of i7(/u) and the relevant 
tropical analogue of the Riemann theta function. We will call them 'tropical . . .' rather casually 
without identifying an underlying tropical geometric objects hoping not to cause a too much 
embarrassment . 

3. In [28], the dynamical system on Cg) • • • (g) i?^''^ equipped with the unique time evolu- 
tion T^^ of the form (I2l6]) was studied under the name of the generalized periodic box-ball 
system. See Sections 12.11 and 12.21 for the notations S^'' and Tj^^^ . Approaches by ultradiscretiza- 
tion of the periodic Toda lattice also capture the but conceivably {T^^^ | / > 1} at most. 
Our periodic An^ SCA is the generalization of the generalized periodic box-ball system with 
V/j = 1 that is furnished with the commuting family of time evolutions {T^'^^ \ i^r<n,l>l}. 
As we will see in the rest of the paper, it is crucial to consider this wider variety of dynamics and 
their whole joint spectrum {E^^^}. They turn out to be the necessary and sufficient conserved 
quantities to formulate the inverse scattering method. Such a usage of the full family {t/''^} 
was firstly proposed in the periodic setting in [25( I26j . In particular in the latter reference, 
the most general periodic An^ SCA on B'^^'^^ (8) • • • B^^'^^ endowed with the dynamics {T^^^} 
was investigated, and the dynamical period and a phase space volume formula were conjectured 
using some heuristic connection with combinatorial Bethe ansatz. This paper concerns the basic 
case Vrfc = VZ^ = 1 only but goes deeper to explore the linearization scheme under which the 
earlier conjectures \25\ [26] get refined and become simple corollaries. We expect that essential 
features in the inverse scattering formalism are not too much influenced by the choice of {r^,, Z^}. 

1.4 Contents of paper 

Let us digest the contents of the paper along Sections [2]-[H 

In Section [21 we formulate the periodic An^ SCA. It is a solvable vertex model [1] associated 
with the quantum affine algebra Ug{An^) at g = in the sense that notions concerning Uq{An^) 
are replaced by those from the crystal theory |14[ [T2] 113] . a theory of quantum group at g = 0. 
The correspondence is shown between the first and the second columns of the following table, 
where we hope the third column may be more friendly to the readers not necessarily familiar 
with the crystal theory. 



^It is a Bethe ansatz folklore that "the solution exists before the model is constructed" according to 
V.V. Bazhanov. 
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Uq{An^) vertex model 


SCA {q = 0) 


combinatorial description 


Kinliov-Kesnetikhm module Wi 
quantum R 
fusion transfer matrix 
conserved quantity 


crystal B ' 
combinatorial R 
time evolution 

(r) 

energy Ej^ 


1 7 "j-llj-ll 

shape r x I semistandard tableaux 
rule (|2.ip on Schensted insertions 
diagram (pr2]) with v = v' £ B""'^ 
n-tuple of Young diagrams /i 



The definition of the Kirillov-Reshetikhin module was firstly given in [191 Definition 1.1], 
although we do not use it in this paper. As mentioned previously, there is the n-tuple of 
commuting series of time evolutions T^'-^^ . . . , T^^"^ with / > 1. The first series t/^^ is the 

ultradiscrete Toda/KP flow [71 [39], and especially its top admits a simple description 
by a ball-moving algorithm like the periodic box-ball system. See Theorem 12.11 and comments 
following it. However, as again said previously, what is more essential in our approach is to make 
use of the entire family and the associated conserved quantities called energy \^E^^^^. 

It is the totality of the joint spectrum {^E^^^^ that makes it possible to characterize the level 
set V{^) ()2.1ip by the n-tuple of Young diagrams and the whole development thereafter. 

As an illustration, Tj^^^ (241123431) = 423141213 is calculated as 

241123431 
12 I 12 I 13 I 11 I 11 I 11 I 12 I 13 I 23 I 12 

34 ~r 23 ~r 24 ~r 24 ~r 24 ~r 22 ~r 23 ~r 24 ~r 34 ~r 34 

4 2 3 1 4 1 2 1 3 

by using the 'hidden variable' called carrier on horizontal edges belonging to i?^'^. This is 
a conventional diagram representing the row transfer matrix of a vertex model [1] whose auxiliary 

(horizontal) space has the 'fusion type' B^'"^. The general case T^^^^ is similarly defined by 
using i?'"''. The peculiarity as a vertex model is that there is no thermal fluctuation due to 
the crystallizing choice q = Q resulting in a deterministic dynamics. Note that the carrier has 
been chosen specially so that the leftmost and the rightmost ones coincide to match the periodic 
boundary condition. This non-local postulate makes the well-definedness of the dynamics highly 
nontrivial and is in fact a source of the most intriguing features of our SCA. It is an important 
problem to characterize the situation in which all the time evolutions act stably. With regard 
to this, we propose a neat sufficient condition in (j2.14p under which we will mainly work. See 
Conjecture 12. 11 Our most general claim is an elaborate one in Conjecture 13.41 In Section [2.31 we 
explain the rigged configurations and Kerov-Kirillov-Reshetikhin bijection which are essential 
tools from the combinatorial Bethe ansatz [15[ I16j at q = 1. Based on Theorem 12.21 we identify 
the solitons and strings in (|2.19p . 

In Section [31 we present the inverse scattering formalism, the solution algorithm of the initial 
value problem together with applications to the volume formula of the phase space and the 
dynamical period. The content is based on a couple of conjectures whose status is summarized 
in Section 13.81 Our strategy is a synthesis of the combinatorial Bethe ansatze at g = 1 [15[ [T6] 
and g = |19j . or in other words, a modification of the KKR theory to match the periodic 
boundary condition. The action-angle variables are constructed from the rigged configurations 
invented at g = 1 by a quasi-periodic extension (j3.ip followed by an identification compatible 
with the q = Q Bethe equation (|3.42p . More concretely, the action variable is the n-tuple of 
Young diagrams ^ = ■ • • , /i^"^) preserved under time evolutions. The angle variables live in 
a tropical analogue of the Jacobi-variety i7(/u) and undergo a straight motion with the velocity 
corresponding to a given time evolution. Roughly speaking, the action and angle variables 
represent the amplitudes and the phases of solitons, respectively as demonstrated in Example ll.il 
The modified KKR bijection $, yield the direct and inverse scattering maps. Schematically 
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these aspects are summarized in the commutative diagram (Conjecture 13. 2p : 



T 



r 



where T stands for the commuting family of time evolutions {T^^''^}. Since its action on 

is linear, this diagram achieves the solution of the initial value problem conceptually. Practical 

calculations can be found in Examples 13.21 and 13. 5[ 

In Section [3T3l we decompose the level set V{fi) further into connected components, i.e., T- 
orbits, and identify each of them with T^^/F-yZ^, the set of integer points on a torus with g 
and F-y explicitly specified in ()2.9p and (j3.25p . Here 7 denotes the order of symmetry in the 
angle variables. As a result we obtain (Theorem 13. ip 



|A („)(mf\pS'^))| . />),^W_1 

\V{,)\ =^ detF, n ia), (.) = (deti^) U (.) 

' V ' 

number of 7"-orbits 

The first equality follows from the decomposition into tori and the second one is a slight calcu- 
lation. The last expression was known as the number of Bethe roots at g = [19]. Thus this 
identity offers the Bethe ansatz formula a most elegant interpretation by the structure of the 
phase space of the periodic SCA. 

Once the linearization scheme is formulated, it is straightforward to determine the dynamical 
period, the smallest positive integer M satisfying T-^ {p) = p for any time evolution T £ T and 
path p G V{fi). The result is given in Theorem 13.21 and Remark 13.41 We emphasize that (13.39P 
is a closed formula that gives the exact (not a multiple of) dynamical period even when there 
are more than one solitons with equal amplitudes and their order of symmetry 7 is nontrivial. 
See Example 13.71 in this paper for n = 2 and also j24| Example 4.10] for n = 1. In Section [3.61 
we explain the precise relation of our linearization scheme to the Bethe ansatz at g = |19] . 
The angle variables are actually in one to one correspondence with what we call the Bethe root 
at g = 0. These results are natural generalizations of the n = 1 case proved in |241 136] . In 
Section 13.71 we treat the general case (I3.46P relaxing the condition ()2.14p . It turns out that the 
linearization scheme remains the same provided one discards some time evolutions and restricts 
the dynamics to a subgroup T' of T- This uncovers a new feature at n > 1. 

In Sectional we derive an explicit formula for the path p S V{^J,) that corresponds to a given 
action-angle variable (Theorem 14. 3p . This is a tropical analogue of the Jacobi inversion problem, 
and the result is indeed expressed by a tropical analogue of the Riemann theta function (j4.14p : 

6(z) = - min | i *n5n + *zn| 

having the quasi-periodicity 6(z + v) = 6(z) + *vi?~^(z -|- ^) for v G BZ'^. Here the G x G 
period matrix B is specified by (j4.6p from the n-tuple of Young diagrams ^ as in (j2.5p - (j2.9p . 
Theorems 14.21 and 14.31 are derived from the explicit formula of the KKR map by the tropical tau 
function Ttrop 

for A^n^ [23|. It is a piecewise linear function on a rigged configuration related to 
its charge. The key to the derivation is the identity 

6= lim (rtrop(R'C) — divergent part) (1-1) 

RC— >oo 

first discovered in [5D] for n = 1. Here the limit sends the rigged configuration RC into the 
infinitely large rigged configuration obtained by the quasi-periodic extension. See (j4.13p for the 
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precise form. It is known [23] that Ttrop is indeed the tropical analogue (ultradiscrete limit) of 
a tau function in the KP hierarchy [llj . Thus our result can be viewed as a fermionization 
of the Bethe ansatz and quasi-periodic solitons at a combinatorial level. As applications, joint 
eigenvectors of {7]^^^} are constructed that possess every aspect as the Bethe eigenvectors at 
(7 = (Section 14. 2p . the dynamical period is linked with the q = Bethe eigenvalue (Section 14. 3p . 
and miscellaneous calculations of some time average are presented (Section 14. 4p . 

The main text is followed by 4 appendices. App endix lAl recalls the row and column insertions 
following [5] which is necessary to understand the rule ()2.ip that governs the local dynamics. 
Appendix [B] contains a proof of Theorem 12.11 based on crystal theory. Appendix [C] is a quickest 
and self-contained exposition of the algorithm for the KKR bijection. Appendix iDl is a proof of 
Theorem 12.21 using a result by Sakamoto [33] on energy of paths. 

2 Periodic A^^ soliton cellular automaton 

2.1 Crystal B^''^ and combinatorial R 

Let n be a positive integer. For any pair of integers r, I satisfying 1 < r < n and Z > 1, we denote 
by B''^ the set of all r X I rectangular semistandard tableaux with entries chosen from the set 
{1, 2, . . . ,n+l}. For example when n = 2, one has B^'^ = {1, 2, 3}, 5^'^ = {11, 12, 13, 22, 23, 33} 
and B^''^ = { 22 ' 23 ' 33 ' 23 ' 33 ' 33 } ■ '^^^ equipped with the structure of crystal [HI [34] 

for the Kirillov-Reshetikhin module of the quantum affine algebra Uq(^An^^- 

There is a canonical bijection, the isomorphism of crystals, between B^'^^B^'^ and B^'^^B^'^ 
called combinatorial R. It is denoted either by R{b c) = d (db' or R{c! ®b') = b® c, or simply 
b ® c d ®b' , where 5 ® c € i?'"'' Cg) B^'^ and d ®b' ^ B^'^ (g) B^'K R is uniquely determined by 
the condition that the product tableaux c • b and b' ■ d coincide [3l]. Here, the product c • 5 for 
example signifies the column insertion of c into 6, which is also obtained by the row insertion 
of b into c [5]. Since c and d are single numbers in our case, it is simplest to demand the equality 



to find the image 6 (8) c ~ d ®b' . See Appendix |A] for the definitions of the row and column 
insertions. The insertion procedure also determines an integer H{b® d) = H[d ® b') called local 
energy. We specify it as -ff = or 1 according as the shape of the common product tableau is 
((/ + l),r-^) or (r,l), respectiveljH. Note that R and H refer to the pair (r, I) although we 
suppress the dependence on it in the notation. 

Example 2.1. Let r = 2, / = 3. 




(2.1) 




111 

223 



product tableau 



1112 
223 



F = 




product tableau 



112 
223 
3 



H = 1. 



We depict the relation R{b c) 



d ®b' (HU and H{b c) = e as 



c 



b 



e 



^Up to a constant shift, this is the original definition of local energy [12] times (— 1). 
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We will often suppress e. The horizontal and vertical lines here carry an element from B^' 
and B^'^, respectively. We remark that R is trivial, namely R{b ® c) = 6 (8) c, for (r, I) = (1, 1) 
by the definition. 



2.2 Definition of dynamics 

Fix a positive integer L and set B = [B^'^)®^. An element of B is called a path. A path 
6i (g) • • • (gi 6i will often be written simply as a word 6162 ■ ■ - bL- Our periodic An'^ soliton cellular 
automaton (SCA) is a dynamical system on a subset of B. To define the time evolution Ti''^ 
associated with B^' , we consider a bijective map B"'^^ (g) B ^ B (g) S*^'^ and the local energy 
obtained by repeated use of R. Schematically, the map 

S*"'' g)B B(g) B^^^ 

■u (g) (61 (g) • • • (g) 6l) (g) • • • (g) h'^) (g v' 

and the local energy ei, . . . , ei, are defined by the composition of the previous diagram as follows. 

hi 62 bi 



ei 


62 




I 


4 I 





■v' 



h'^ (2-2) 

The element v or v' is called a carrier. Set = 6'^(g - • - (gft^. Given a path p = hi<gi- ■ -igbL & B, 
we regard v' = v'{v;p),ek = ek{v;p) and p' = p'{v;p) as the functions of v containing p as 
a 'parameter'. Naively, we wish to define the time evolution T^^^ of the path p as Ti^\p) = p' 
by using a carrier v that satisfies the periodic boundary condition v = v'{v;p). This idea indeed 
works without a difficulty for n = 1 |24pl . In addition, T^^^ is always well defined for general n, 
yielding the cyclic shift: 

ri^Hbi ig • • • (g bi-i (g 6l) = 6i (g 61 (g • • • (g hL-i- (2.3) 

This is due to the triviality of R on B^'^ (g mentioned above. In fact the unique carrier is 
specified as v = v' = b^. Apart from this however, one encounters the three problems 
to overcome in general for n > 2. For some p £ B, one may suffer from 

(i) Non-existence. There may be no carrier v satisfying v = v'{v;p). 

[ii) Non- uniqueness. There may be more than one carriers, say vi,V2,... ,Vm such that Vj = 
v'{vj;p) but p'{vi\p) ^p'{vj;p) or ek{vi]p) / ek{vj;p) for some i, j, k. 

To cope with these problems, we introduce the notion of evolvability of a path according to |26j . 
A path p £ B is said Tj^^^-evolvable if there exist v G B^'^ such that v = v'{v;p), and moreover 
p' = p'{v;p), Ck = ek{v;p) are unique for possibly non-unique choice of In this case we define 

t['\p)=p', E\'\p) = ei + --- + eL, (2.4) 

indicating that the time evolution operator is acting on p. The quantity Ej^ [p) € Z>o is 
called an energy of the path p. 

By the standard argument using the Yang-Baxter equation of the combinatorial R, one can 
show 

*Even for n = 1, it becomes nontrivial for higher spin case |22j . 

^We expect that the uniqueness of the local energy ek{v;p) follows from the uniqueness of p'{v;p) alone. 
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Proposition 2.1 (|26j). The commutativity T^"''^T^^\p) = Tj^^T-"'\p) is valid if all the time 



evolutions T^''^ here are acting on T^^^ -evolvahle paths. Moreover, E^""^ (T^''\p)) = e\"'\p) and 



Ef\T!>''\p)) = Ef\p) hold. 

(r) 

If p is -evolvable for all r, /, then it is simply called evolvable. The third problem in 
defining the dynamics is 

[Hi) Even if p is evolvable, T^^^p) is not necessarily so in general. 

For instance, p = 112233 G {B^^^)®^ is evolvable but ^(p) = 213213 is not ^-evolvable. 
In fact, the non-uniqueness problem [ii) takes place as follows: 

213213 213213 
11122112211112 



32233333333223 
311223 223311 



To discuss the issue {Hi), we need to prepare several definitions describing the spectrum ^^E^^. 
Let /i= be an ?i-tuple of Young diagrams. From we specify the data 



(a) 



9a as in the left diagram in (j2.5p . 



(a) 



(a) 



1^ ha I 



m- 



(a) 



1^ ha I 



(a) 

f 1 

1,1 



Here, m\^'^ x if^'^ rectangle part is called the (a, i) block. Furthermore we introduc^ 

H = \^{a,i,a) \ 1 < a < n, I < i < ga, 1 < a < 
H = {{a,i) |l<a<n, l<i< ga}, 

pf^=L6a,- Ca,mm{ll^\lf^)=L6ai- Yl Ca,ui\n{lf ,lf)mf , 

{bjl3)(iH {bj)eH 

F = {Fai,bj)(^aiUbj)&H^ Fat,bj = ^ab^tjpf^ + Cab min {if ,lf)mf , 

n Qa n 



(2.5) 

(2.6) 
(2.7) 

(2.8) 
(2.9) 



a=l i=l 



where (C^ 



) l<a,b<n 



is the Cartan matrix of An, i.e.. Cab = "^^ab — Sa,b+i — Sa,b-i- The quantity 



p\ is called a vacancy number. 

®The meaning of i in m'f'' and p*°' here has been altered(!) from the hteratures [l9l gOl [211 [22l [23l IM [251126] . 
When the space is tight, we will often write (aia) for (a, i,a) £ H and (ai) for (a, i) £ H, etc. 
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Given ^ = {fi^ define the sets of paths B D V D V{fj,) D Vj^[^) by 

V = {p^B\ #(1) > #(2) > • • • > #(n + 1)}, (2.10) 

V{p) = |p G P I p : evolvable, sj"^(p) = ^ min (2.11) 

V +{^l) = {peV ill) \ p:\iig\iesi}, (2.12) 

where #(a) is the number of a G B^'^ inp = bi® ■ ■ ■ ® bi- In view of the Weyl group symmetry 
|261 Theorem 2.2], and the obvious property that T^^^'^ is weight preserving, we restrict ourselves 
to V which is the set of paths with nonnegative weight. A path hi ® ■ ■ ■ ® bL is highest if the 
prefix 6i (g) • • • (g) 6j satisfies the condition #(1) > #(2) > • • • > #(n + 1) for all 1 < j < L. We 
call V{fi) the level set associated with /x. The n-tuple /i of Young diagrams or equivalently the 
data ("T'i'^\ 'i°^)(aj)g7f will be referred to as the soliton content of V{fi) or the paths contained 
in it. This terminology comes from the fact that when L is large and #(1) ^ ^(2), . . . , ^(n + 1), 
it is known that the paths in P(//) consist of 7n^^^ solitons with amplitude separated from 
each other or in the course of collisions. Here a soliton with amplitude I is a part of a path of 
the form ji (g) j2 (E> • • • 'g' ji G (iJ^'^)®' with ji >■■■> ji > 2 hy regarding 1 G B^'^ as an empty 
background. Note that a soliton is endowed not only with the amplitude / but also the internal 
degrees of freedom {ji}. The data m[°'\ l^f^ with a > 1 are relevant to the internal degrees of 
freedom of solitons. 

Note that the relation in ()2.1ip can be inverted as {E^^ := 0) 

-E^\ + 2e\'^^ - Eil\ = h'^'' ^ = ^ 1 < ^ < (2.13) 

I otherwise. 

Therefore the correspondence between the soliton content = (/x^^^ . . . ,^^"^) and the energy 

{E^°^} is one to one. Either or /i are conserved quantities in the following sense; if p G 

V{ijl) and T^^'\p) is evolvable, then Proposition 12.11 tells that Tj^'^\p) G Vip). Comparing (|2.7p 
and (|2.1ip . one can express the vacancy number also as 

pf^=L5^i-jZCa,E%. 
b=i 

Now we are ready to propose a sufficient condition under which the annoying feature in 
problem {Hi) is absent: 

= (^(^\ . . . satisfies pf > 1 for ah {a,i) G H. (2.14) 

This condition was first introduced in |22j. 

Conjecture 2.1. Under the condition (j2.14p . T^^\v{fi)) = V{fi) holds for any r, I. 
For a most general claim, see Conjecture [32 

for arbitrary times on any path p G 'P{^). In particular the inverse (T)^^^) ^ exists. We denote 
by S(p) the set of all paths generated from p in such a manner. By the definition it is a subset 
of the level set V{^). Call S(p) the connected component of the level set containing the path p. 

One can give another characterization of S(p) in terms of group actions. Let T be the 
abelian group generated by all T^^^^s. Then T acts on the level set Vifi)- This group action 
is not transitive in general, i.e., V{fi) is generically decomposed into several T-orbits. In what 
follows, connected components and T-orbits mean the same thing. 



The equality T^^\v{ii)) = 'P(/i) implies that one can apply any time evolution in any order 
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Example 2.2. Let n = 2,L = 8 and /i = ((211), (1)). The vacancy numbers are given as 
p"^^ = 1, p"2^ = 3, p^^^ = 1, hence the condition ()2.14p is satisfied. 

m □ 



Then we have P(/i) = T.{px) U S(py) with px = 11221123 and py = 11221213. The connected 
components are written as S(px) = |J^=i -^i and S(py) = |Ji*=i where Xj and Yi are given in 
the following table. 



i 


Xi Yi 


1 


[11221123] 


[11221213] 


2 


[11122123] 


[11121223] 


3 


[21112123] 


[21121123] 


4 


[21121213] 


[21211213] 


5 


[21212113] 


[12212113] 


6 


[12122113] 


[11212213] 


7 


[12112213] 


[11211223] 


8 


[12111223] 


[21211123] 


9 


[12211123] 


[12211213] 



Here the symbol [ ] denotes the set of all cyclic shifts of the entry. For instance, [123] = 
{123,231,312}. For any i G Z/9Z we have Ti^\Xi) = Xi, T^^^iXi) = Xi+i, T^I{Xi) = Xi+3 
and the same relations for Yi. Hence the claim of Conjecture 12.11 is valid for this V{p). 

The time evolution T^'^ has an especially simple description, which we shall now explain. 
Let Bi be the set of paths in which 1 G B^'^ is contained most. Namely, 

Bi = {peB\ tt(l) > S(a), 2 < a < 71 + 1}. (2.15) 

Thus we see V C Bi C B. We define the weight preserving map Ka : i?i — )• -Bi for a = 
2, 3, . . . , n + 1. Given a path p S i?i regarded as a word p = 6162 • • • &l G {1, . . . , n + 1}-^, the 
image Ka{p) is determined by the following procedure. 

(i) Ignore all the numbers except 1 and a. 

(ii) Connect every adjacent pair al (not la) by an arc, where 1 is on the right of a cyclically, 
[in) Repeat [ii) ignoring the already connected pairs until all a's are connected to some 1. 
[iv) Exchange a and 1 within each connected pair. 

Theorem 2.1. Suppose p £ Bi and p is T^^ -evolvahle. Then, is factorized as follows: 

T^^\p) = K2K3---Kn+i{p). (2.16) 

A proof is available in Appendix [Bl where each Ka is obtained as a gauge transformed simple 
reflection for type A^^ affine Weyl group. The dynamics of the form (j2.16p has an interpretation 
as the ultradiscrete Toda or KP flow [3 [39]. See also [281 Theorem VI. 2]. As an additional 
remark, a factorized time evolution similar to ()2.16p or equivalently (IB.4P has been formulated 
for a periodic SCA associated to any non-exceptional affine Lie algebras [25]. In the infinite 
(non-periodic) lattice case, it was first invented by Takahashi [37] for type The origin of 

such a factorization is the factorization of the combinatorial R itself in a certain asymptotic 
domain, which has been proved uniformly for all non-exceptional types in [8]. 
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Example 2.3. In the left case of Example ll.il the time evolution Tg^"* is actually equal to 

.(1) _ ^(1) 



The time evolution from i = 2 to t = 3 is reproduced by To 



K2K^ as follows: 



n I r~n I — I r 

331132111321111221122213 



1 



n I M n I 
113312331123111221122211 



213311332113211112211122 



K2 



The procedure {i)—{iv) is delightfully simple but inevitably non-local as the result of expelling 
the carrier out from the description. It is an interesting question whether the other typical time 
evolutions , . . . , T^'^ admit a similar description without a carrier. 

2.3 Rigged configuration 

An n-tuple of Young diagrams /i = ■ ■ ■ , /i^"^) satisfying pf^^ > for all (a, i) £ H is called 
a configuration. Thus, those fi satisfying (I2.14p form a subset of configurations. Consider the 
configuration attached with the integer arrays called rigging r = (r("))i<a<„ = {rl"']^)i^aia)&H 
in the right diagram in (12.50 . The combined data ()U,r) = {{fi^^\r^^^), . . . , {fi^'^\r^^^)) is called 
a rigged configuration if the condition 

< rjf < < • • • < r^"^ < ^ for aU (a, i) eH (2.17) 

is satisfieclll. We let RC(/x) denote the set of rigged configurations whose configuration is fi. It 
is well known \15\ [16] that the highest paths in B are in one to one correspondence with the 
rigged configurations by the Kerov-Kirillov-Reshetikhin (KKR) map 

(l):{peB\p: highest} ^ U^RC(^), (2.18) 

where the union is taken over all the configurations. The both cj) and (/)~^ can be described 
by an explicit algorithm as described in Appendix O Our main claim in this subsection is the 
following, which is a refinement of ()2.18p with respect to // and an adaptation to the periodic 
setting. 

Theorem 2.2. The restriction of the KKR map (p to V+{fi) separates the image according to fx: 

The proof is available in Appendix[Dj We expect that this restricted injection is still a bijec- 
tion but we do not need this fact in this paper. 

Let us explain some background and significance of Theorem 12.21 Rigged configurations 
were invented [H] as combinatorial substitutes of solutions to the Bethe equation under string 
hypothesis [2]. The KKR map is the combinatorial analogue of producing the Bethe vector 
from the solutions to the Bethe equation. The relevant integrable system is sln+i Heisenberg 
chain, or more generally the rational vertex model associated with Uq[An^^ at g = 1. In this 
context, the configuration /i (j2.5p is the string content specifying that there are in^'^'^ strings 

^We do not include L in the definition of rigged configuration understanding that it is fixed. It is actuaUy 
necessary to determine the vacancy number p'"' (|2.7|l appearing as the upper bound of the rigging (|2.17|) . 
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with color a and length l^^ . Theorem 12.21 identifies the two meanings of Namely, the soliton 
content for V{^) measured by energy is equal to the string content for RC(/x) determined by the 
KKR bijection. Symbolically, we have the identity 

soliton = string, (2.19) 

which lies at the heart of the whole combinatorial Bethe ansatz approaches to the soliton cellular 
automata [20l[2ll[22l[23l[2ll[25l[26]. It connects the energy and the configuration by (|2T3]1 . 
and in a broader sense, crystal theory and Bethe ansatz. Further arguments will be given 
around (j4.23p . For n = 1, Theorem 12.21 has been obtained in \14\ Proposition 3.4]. In the sequel, 
we shall call the n-tuple of Young diagrams /i = {^jf'^\ ■ ■ ■ , n^"'^) either as soliton content, string 
content or configuration when p^""^ > for all (a, i) E H. 

A consequence of Theorem 12.21 is that for any path bi® ■ ■ ■ G V{^) with soliton content 

the number of a G B^'^ in bi, ... ,bi is given by 

#(a) = - (l<a<?i + l), (2.20) 

where we set l^*-"-*! = L and j^'-""''^^! = 0. This is due to a known property of the KKR map (p 
and the fact that fi is a conserved quantity and Conjecture 13.51 

Example 2.4. Let n = 2,L = 8 and consider the same configuration fi = ((211), (1)) as 
Example 12.21 The riggings are to obey the conditions < r^il,r^i < 1 and < r^l < r^l < 3. 

(/i(2),r(2)) 



Hence there are 2 • 2 • 10 = 40 rigged configurations in RC(/u). The elements of 'P+(/x) and the 
riggings for their images under the KKR map (p are given in the following table. 



i 


highest paths in Xi 


riggings 


highest paths in Yi 


riggings 


1 


11221123, 11231122, 12311221 


0331, 1110, 0000 


11221213, 12131122 


0321, 1100 


2 


11122123, 11221231, 12311122 


1331, 0221, 1000 


11121223, 11212231 


1321, 0211 


3 


11121232, 11212321 


1221, 0111 


11211232, 12112321, 11232112 


1211, 0101, 0310 


4 


11212132, 12121321, 12132112 


nil, 0001, 0300 


12112132, 11213212, 12132121 


1101, 1310, 0200 


5 


12121132, 12113212 


1001, 1300 


12113122 


1200 


6 








11212213 


0311 


7 


12112213, 11221312 


0301, 0330 


11211223, 12112231, 11223112 


1311, 0201, 0320 


8 


12111223, 11122312, 11223121 


1301, 1330, 0220 


12111232, 11123212, 11232121 


1201, 1320, 0210 


9 


11123122, 11231221 


1220, 0110 


12131221, 11213122 


0100, 1210 



Here the riggings denote r\ Ir 2 



(/iW,r«) 



P^' 



.(1) 
1,1 



,(1) 
2,2 

,(1) 
2,1 



3 Inverse scattering method 
3.1 Action-angle variables 

By action-angle variables for the periodic An^ SCA, we mean the variables or combinatorial ob- 
jects that are conserved (action) or growing linearly (angle) under the commuting family of time 
evolutions {t/''^}. They are scattering data in the context of inverse scattering method [6l [3l H] . 
In our approach, the action-angle variables are constructed by a suitable extension of rigged 
configurations, which exploits a connection to the combinatorial Bethe ansatze both at g = 1 [15] 
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and q = [19]. These features have been fully worked out in |24j for More recently it has 
been shown further that the set of angle variables can be decomposed into connected compo- 
nents and every such component is a torus [36]. Here we present a conjectural generalization of 
these results to An^ case. It provides a conceptual explanation of the dynamical period and the 
state counting formula proposed in [25} [26]. 

Consider the level set V{^) (12. IIP with fx = {^^^\ . . . satisfying the condition ()2.14p . 

The action variable for any path p G V{fJ,) is defined to be itself. By Proposition 12. 1|, it is 
a conserved quantity under any time evolution. 

Recall that each block (a, i) of a rigged configuration is assigned with the rigging r^-"^ with 

a = 1, . . . , m\"'\ We extend r^^*^ to a G Z uniquely so that the quasi-periodicity 

r^^^^=rS.+I^ iaeZ) (3.1) 
is fulfilled. Such a sequence r = (rj^)^^^ will be called the quasi-periodic extension of 

Set 

{ai)eH 

A{m,p) = {{Xa)aez I A„ G Z, Aq, < K+i, Xa+m = K + P for all a} (m > l,p > 0), (3.2) 
where Y\ denotes a direct product of sets. Using Theorem 12. 2[ we define an injection 

.:P+(/.) ArC(/.) Jifi) 

p+ (M,r) ^ (^•a)(ai)eF,aeZ' 

where the bottom right is the quasi-periodic extension of the original rigging r = (^i q) (^j^^jg^:^ 

in (fi, r) as explained above. Elements of i7(/i) will be called extended rigged configurations. 
One may still view an extended rigged configuration as the right diagram in ()2.5p . The only 
difference from the original rigged configuration is that the riggings rj^ outside the range 
1 < a < nrif^ have also been fixed by the quasi-periodicity (j3.ip . In what follows, either the 
rigging or its quasi-periodic extension („j)g;H',aGZ ^ '^^^^ denoted by 

the same symbol r. 

Let T be the abelian group generated by all the time evolutions T^"''^ with 1 < a < n and 

I > 1. Let further ^ be a free abelian group generated by the symbols sj"^ with (a,i) G H. We 
consider their commutative actions on i7(/x) as follows: 



^l'^^ ■■ i^tl)(M)eH.BeZ ^ {''fi + ^-b rain {I, if (3.3) 



■ i^fi)m.H,^.z ^ ^!l-^S..S, + Ga^nrin ^ ))(,,),^,,,,. (3.4) 
The generator s-"^ is the An^ version of the slide introduced in [23] for n = 1. Now we define 



Jil^) = Jifi)/A, (3.5) 

which is the set of all ^-orbits, whose elements are written as ^ • r with r G Elements of 

will be called angle variables. See also (j4.17p . Since T and A act on commutatively, 
there is a natural action of T on J{n). For t £ T and y = A • r £ JifJ-), the action is given by 
t ■ y = A ■ {t ■ r). In short, the time evolution of extended rigged configurations (|3.3p naturally 
induces the time evolution of angle variables. 
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Remark 3.1. From (j3.4p and ()2.7p . one can easily check 




for r = {ff^p) {pj-^^jj i^^i £ Therefore r and (^^^^^ {v) define the same angle variable. 

3.2 Linearization of time evolution 

Given the level set V{^) with ^ satisfying the condition (j2.14p . we are going to introduce the 
bijection <I> to the set of angle variables Recall that denotes the connected component 

of P(/u) that involves the path p. 

Conjecture 3.1. Under the condition (j'2.14p . n'P+(/.f) 7^ holds for any p E V{p). 

This implies that any path p can be expressed in the form p = t • by some highest path 
p+ G V+{fi) and time evolution t = JK^/''^)'^' ^ 

Although such an expression is not unique 
in general, one can go from 'P(^) to along the following scheme: 

P(/i) TxV+if,) ^ rx J(^) J(^) ^gg^ 

p ^ {t,p+) ^ (t,r) ^^.(t-r). 

Note that the abelian group T of time evolutions acts on paths by (|2.4p and on extended rigged 
configurations by ()3.3p . The composition (j3.6p serves as a definition of a map $ only if the 
non-uniqueness of the decomposition p 1— )• is 'canceled' by regarding t ■ r mod ^. 

Our main conjecture in this subsection is the following. 

Conjecture 3.2. Suppose the condition ()2.14p is satisfied. Then the map ^ ()3.6p is well-defined, 
bijective and commutative with the action of T. Namely, the following diagram is commutative 



T 



r (3.7) 



Tip.) Jill) 

This is consistent with the obvious periodicity {T^^)^{p) = p of paths under the cyclic 
shift (12. 3p and Remark 13.11 The time evolution of the angle variable (13. 3p is linear. Thus the 
commutative diagram (|3.7p transforms the nonlinear dynamics of paths into straight motions 
with various velocity vectors. The set J{fi) is a tropical analogue of Jacobi variety in the theory 
of quasi-periodic solutions of soliton equations [31 H]. The map <I> whose essential ingredient is 
the KKR bijection (p plays the role of the Abel-Jacobi map. Conceptually, the solution of the 
initial value problem is simply stated as t^{p) = o t-^ o for any t £ T, where t in 

the right hand side is linear. Practical calculations can be found in Examples 13.21 and 13.51 In 
Section [H we shall present an explicit formula for the tropical Jacob i-inversion in terms 
of a tropical Riemann theta function. The result of this sort was first formulated and proved 
in pm^EU forn = 1. 

Example 3.1. Take an evolvable path 

p = 211332111321133112221112 G (S^'^)®^^, (3.8) 



^We do not restrict f to a specific subgroup of T. 
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whose energies are given as 



£;W = 5, 4^^ = 10, 4^j = 12, e['^=2, =3, i^f =4, Sg = 5. 

From this and we find p G with = ((33222), (41)). The cyclic shift {T^^^y{p) 

is not highest for any j. However it can be made into a highest path and transformed to 
a rigged configuration as follows: 



(/i(i),rW) 



(/i(2),r(2)) 



= 111221113221132113311322 



















6 
5 
1 











4 2 
2 1 














(3.9) 



Thus action variable (/i^^^/x^^^) indeed coincides with fi = ((33222), (41)) obtained from the 
energy in agreement with Theorem 12.21 Although it is not necessary, we have exhibited the 
vacancy numbers on the left of fi^^^ and /i^^^ for convenience. We list the relevant data: 



(ai) 




(a) 

m- 


a) 

Pi 




(a) 

Ti 


(11) 


3 


2 


4 


2,4 


2 


(12) 


2 


3 


7 


1,5,6 


1 


(21) 


4 


1 


2 





1 


(22) 


1 


1 


1 





1 



The index set H (12. 6p corresponding io p. = fi^"^^) in ()3.9p is 
H = {(111), (112), (121), (122), (123), (211), (221)}, 



(3.10) 



which has the cardinality G = 7. The quantity 7^. ' is the order of symmetry which will be 
explained in Section [3.31 The matrix F ()2.8p reads 



(1) 



4mi^^ 



4m2^^ 



(2) 



-2m, 



-m 



-m, 



W 
2 

(1) 



-2mf^ 



— m. 
— m, 
2m 



2m 



(2) ^(2) 



(2) 
2 

(2) 
2 

(2) 
2 

(2) 



/16 


12 


-3 


-1\ 


8 


19 


-2 


-1 


-6 


-6 


10 


2 


\-2 


-3 


2 


3/ 



(3.11) 



We understand that the riggings here are parts of extended one obeying (13. ip . Then from (13.91 
and ()3.3p . the extended rigged configuration for p ()3.8p is the left hand side of 



1 



,(1) 



















9 
7 

6 











(3.12) 
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Here the transformation to the right hand side demonstrates an example of identification by A. 
These two objects are examples of representative elements of the angle variable for p 



We wiU calculate (t/^V^^Ip) in Example [33] using the reduced angle variable that will be 
introduced in the next subsection. 

Example 3.2. Consider the evolvable path 

321113211222111223331111 G (bI.I)®^'' 



P 



(3.13) 



which is same as in Example 11.11 The cyclic shift {T^^^y{p) is not highest for any j. However 
it can be made into a highest path and transformed to a rigged configuration as 



p+ = (r«)2irf)(p) 

= 111221322111221332111331 



















6 
3 
1 











1 


2 













1 








(3.14) 



Thus p (I3.13P belongs to the same level set Vin) with /x = ((33222), (41)) as Example 13. H 
hence the matrix F (j2.8p is again given by (j3.1ip . However, it has a different (trivial) order of 
symmetry 7^"^ given in the following table: 



(ai) 




(a) 

m- 


(a 
Pi 




(a) 


(11) 


3 


2 


4 


0,1 


1 


(12) 


2 


3 


7 


1,3,6 


1 


(21) 


4 


1 


2 


2 


1 


(22) 


1 


1 


1 





1 



Regard (j3.14p as a part of extended rigged configuration as in Example 13. 1[ Then from (j3.3p , 
the extended rigged configuration for p (I3.13P is the left hand side of 



20 










21 






1 



,(2) 



-15 
-18 
-20 



23 










24 




1 



11 



-17 
-20 
-22 



(3.15) 



Here again the transformation to the right hand side demonstrates an example of identification 
by A. These two objects are examples of representative elements of the angle variable for p (j3.13p . 
To illustrate the solution of the initial value problem along the inverse scheme (|3.7p , we derive 



221132211331111321322111. 



(3.16) 



By applying (13. 3p to the left hand side of (I3.15p . the angle variable for (T, 
as 



ji)^iooo^^-j is obtained 
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1985 
1982 
1980 
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where the equivalence is seen by using (s^ 

hand side is the angle variable of (T^^Y T'^\p'_y_) , where p'_^ is the highest path obtained by the 
KKR map 
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112211132112221331133211 



Calculating (T^^^^) one finds (IXTHD . 

3.3 Decomposition into connected components 

As remarked under Conjecture 12.11 the level set V{^) is a disjoint union of several connected 
components under the time evolution 7~. It is natural to decompose the commutative 

diagram (j3.7p further into those connected components, i.e., T-orbits, and seek the counterpart 
of in To do this one needs a precise description of the internal symmetry of angle 

variables reflecting a certain commensurability of soliton configuration in a path. In addition, it 
is necessary to separate the angle variables into two parts, one recording the internal symmetry 
(denoted by aJhI and the other accounting for the straight motions (denoted by uj). We call 
the latter part reduced angle variables. The point is that action of T becomes transitive by 
switching to the reduced angle variables from angle variables thereby allowing us to describe the 
connected components and their multiplicity explicitly. These results are generalizations of the 
n = 1 case [36] . 

Let m > 1 and p > 00 We introduce the following set by imposing just a simple extra 
condition on A(m,p) (13. 2p : 



A{m,p) = {{Xa)a£Z I Aq G Z, Ai = 0, A„ < A«+i, X^+m = Xa+P for all a}. 
This is a finite set with cardinality 

\A{m,p)\ 

Given a configuration fi, specify the data mf^\ pf \ etc by (j2.5p - (|2.9p . Set 



p + m — 1 
m — 1 



(3.17) 



(3.18) 



\{_K{mf\pt^), 

{ai)&H 

and consider the splitting of extended rigged configurations into and X'^ as follows: 
J{fi) ^ X^ xX'^ 



i > 



(u;,A) 



where uj 



{ai)eH,aeZ 

')f •^^TJ and A = (A^-^^), .^^^7 



i,a i,l ' 



(3.19) 



Recall that T and A are the abelian groups acting on J{fJ.) as in (j3.3p -(f3T 
one correspondence (j3.19p . they also act on X commutatively as 



r/'^)(u;,A) = (a; + /.r,A), 



(a) 



h 



(a) 



(bi)GH' 



By the one to 
(3.20) 



This A may be viewed as an additional conserved quantity. 
^This p is not a path of course and will set to be a vacancy number shortly. 
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In particular, T acts on transitively leaving unchanged. 

From (j3.5p and p.l9p . the angle variables may also be viewed as elements of X/A. It is the 
set of all ^-orbits, whose elements are written as ^ ■ x with x £ X. Since T and A act on X 
commutatively, there is a natural action of T on X/A. For t £ T and y = A-x S X/A, the action 
is given hy t-y = A - {t-x). For any representative element (cj, A) G X/A of an angle variable, we 
call the X^ part uj a reduced angle variable. Under the simple bijective correspondence (j3.19p . 
the map (13. 6p is rephrased as 



P(/i) ^ TxV+ifi) TxX ^ 



p ^ ^ (t,(w,A)) ^ A- {t-{io,\)), 

where we have used the same symbol The conjectural commutative diagram (j3.7p becomes 



r 



r (3.22) 



Now we are ready to describe the internal symmetry of angle variables and the resulting 
decomposition. We introduce a refinement of A(m,p) (j3.17p as follows: 



A^{m,p) = <! A G A ( -, ^ ) A ^ A ( ^, ^ ) for any 7' > 7 |> , (3.23) 



7 7 



m p 



where 7 is a (not necessarily largest) common divisor of m and p. In other words A-y(m,p) is the 
set of all arrays {Xa)a&'L G A{m,p) that satisfy the reduced quasi-periodicity Xa+m/-y = Xa+p/^ 
but do not satisfy the same relation when 7 is replaced by a larger 7'. Such 7 is called the order 
of symmetry. By the definition one has the disjoint union decomposition: 

A{m,p) = y A^(m,p), 

7 

where 7 extends over all the common divisors of m and p. Taking the cardinality of this relation 
using (|3.18p amounts to the identity 



Let 7 = (7i'^^)(aj)g77 be the array of order of symmetry for all blocks. Thus 7^^°^ is any 

nmon divisor 
and X-y C X by 



common divisor of m^""^ and p^""^ . We introduce the g x g matrix and the subsets X^ C X 



F^ = (i^a..,/7f)(.,),(,,),H' (3-25) 

^7 = n {'^t\pt^) c ^^ X.^ = X^ X X2 c X, 

{ai)e7? 

where Fai^bj is specified in (|2.8p . The matrix F hence F^ are positive definite if pf^^ > for any 
(o, i) G H. 
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Example 3.3. Consider Example 13.11 Take the left hand side of (j3.12p as a representative 
element of the angle variable of the path (|3.8p . Then the splitting (|3.19p gives the reduced angle 
variable to = *( — 1, —2, —3, —3) and A is the quasi-periodic extension of the rigging in 



















5 
4 












2 



















(3.26) 



Thus A = {X^i\x^2\ \ A^^^) reads 

aS'^ = (...,-4,-2,0,2,4,6,8, 10, 
Af) = (...,-2,0,2,4,...), 



^2 



A 



with A'"'' = 0. There is order 2 symmetry 7^^^ = 2 in the block (a, i) = (1,1) since p\ 



(2) 



(■ 



-7,-3,-2,0,4,5,7,11,12 
-1,0,1,2,...), 



5 • • • / ' 



(1) 



and m{^^ = 2 have a common divisor 2, and x[^^ satisfies A^^^ = aI^I + p[^^ /2. Thus 



from (|3.1ip . the matrix Fj (j3.25p becomes 
F^ = F-diag(i, 1,1,1) 



/ 8 12-3 -1\ 

4 19 -2 -1 

-3 -6 10 2 

V-1 -3 2 3 / 



(3.27) 



Example 3.4. Consider Example 13.21 Take the left hand side of (j3.15p as a representative 
element of the angle variable of the path (|3.13p . Then the splitting (j3.19p gives the reduced 
angle variable w = *(— 21,— 20,1,— 1) and A is the quasi-periodic extension of the rigging in 



















5 
2 












1 



















Thus A = {X'i\x^^\xf\xf'>) reads 



A^') = (..., -4, -3, 0, 1,4, 5, 8, 9, . . .), A^^^ = (..., -7, -5, -2, 0, 2, 5, 7, 9, 12, . . .), 



(1) 



Af^ = (...,-2,0,2,4,...), 



A, 



(2) 



(...,-1,0,1,2,...) 



with A-"^ = 0. The order of symmetry is trivial in that V7|"^ = 1. In this case one has F~f = F 

dan]). ' 

An important property of the subset X-f C X is that it is still invariant under the actions of 
both T and A. Let X~^/A be the set of all ^-orbits. According to X = jj^-'^^, we have the 
disjoint union decomposition: 



X/A = [_\{X^/A), 



(3.28) 
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which induces the disjoint union decomposition of the level set V{fi) according to (j3.22p . Writing 
the pre-image of X^/A as 'P-y(^), one has V{fJ.) = U^V-f{ij) and the conjectural commutative 
diagram (j3.22p splits into 



T 



T 



(3.29) 



-> X^/A 



V'y{fi) is the subset of the level set V{fi) characterized by the order of symmetry 7. We are 
yet to decompose it or equivalently X^/A further into T-orbits. Note that one can think of an 
^-orbit A ■ X either as an element of Xj/A or as a subset of X~^. Similarly a T-orbit T ■ {A - x) 
in Xj/A can either be regarded as an element of {X~^/A)/T or as a subset of X~^/A. We adopt 
the latter interpretation. Then as we will see shortly in Proposition 13. H the T-orbits T ■ {A - x) 
of X~^/A can be described explicitly in terms of TL^ jF^ Z^, the set of integer points on the torus 
equipped with the following action of T: 



mod 



for I G Z^/F^Z^, 



(3.30) 



where is specified in (j3.20p . In what follows, we shall refer to Z^/F^ Z^ simply as torus. 



Given any x = (w, A) G X~^, consider the map 



X : T-{A-x) 
t-{A-x) 



l^jF^Z^ 



t • ijj 



mod F^Z^. 



(3.31) 



Proposition 3.1 ([36j). The map x is well-defined, bijective and the following diagram is com- 
mutative: 



T- (A-x) 
r 

T- {A-x) 



^ Z9/F^Z3 



T 



7^7 



Combining Proposition 13.11 with Conjecture 13.21 or its refined form (|3.29p . we obtain an 
explicit description of each connected component (T-orbit) as a torus. 

Conjecture 3.3. For any path p G V-y{fi), the map := gives a bijection between the 

connected component S(p) and the torus Z^/F-^Z^ making the following diagram commutative: 



Z9/F^Zy 
T 

J/F^Z9 



(3.32) 



m - 

T 

m - 

The actions of T in Proposition 13.11 and Conjecture 13.31 are both transitive. Conjecture 
is a principal claim in this paper. It says that the reduced angle variables live in the torus 
Z^/F-yZ^, where time evolutions of paths become straight motions. See Conjecture 13.61 for an 
analogous claim in a more general case than ()2.14p . 



Remark 3.2. Due to the transitivity of T-action, the commutative diagram ()3.32p persists even 
if is redefined as + c with any constant vector c £ Z^, meaning that the choice of the 
inverse image of G Z^/F^Z^ is at one's disposal. With this option in mind, an explicit way 
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to construct the bijection is as follows. Fix an arbitrary highest path p+ in and set 
= 0. For any p' € let t be an element of T such that p' = t ■ p^. One can always 

find such t since T acts on T,{p) transitively. If t is written as t = Yiri (^Z^'*)'^' > then we set 

Remark 3.3. For {a,i) G we define r/-°^ = l^*^^ + 1 (i < ga), = 1 (i = ^a). Let T be the free 
abelian group generated by all T^?\ 's. We let it act on the torus Z^/F-y by (13.301) . This action 

is transitive because the (^-dimensional lattice ©(ai)el7 ^^^w C coincides with itself. If 
Conjecture 13.31 is valid, then as a subgroup of T, T also acts on transitively. Hence the 

(r) 

image of the bijection $^ in the previous remark can be written as a linear combination of /i 



with (r, i) & H only. 



Example 3.5. As an application of the inverse scheme (j3.32p . we take the path p ()3.8p in 
Example 13.11 and illustrate the solution of the initial value problem to derive 

(r2^^y°°°(p) = 122211122113321113211331, 
(rj^y°°°(p) = 213321112211112221331113, 
(rf = 331113211332112211111222, 

(r2^^y°°°(p) = 311132113211122211221133, 
(rf = 113221132111221132213311, 

(r]^y°°°(p) = 221113221112211321333111. 

In Example 13. 3|, we have seen that the reduced angle variable of p is given by 



(3.33) 



OJ 



-2 
-3 
V-3/ 



4 

7 ' 



where F-y is specified in (|3.27p . The velocity vector h^""^ (j3.20p of the reduced angle variable 
corresponding to the time evolution t!'"'^ is given by 



.(2) 











^3\ 






1 






2 



, — 




1 










vo/ 




vv 


/o\ 








/ON 












) "'4 — 







2 


3 


4 




vv 




vv 









The higher time evolutions /ij^g and /i^^''^ coincide with /ig^' and /i4^\ respectively. Moreover from 

^j'U)'j2000^^-j^ and the latter is easily obtained since t[ 



one has {t!^^)^^^^ {p) 



is just a cyclic shift ([O]), yielding {T^ ') {p) = 122211122113321113211331 in agreement 
with ([333]) . 
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Next we consider (T.^^^^^^'^^ (p). The reduced angle variable for this path is 



oj + lOOO/i 



(1) 



/2999\ 
1998 
-3 
V-3/ 



/A 

6 

1 

V/ 



mod f = ^ + /if ^ + 



/1\ 







mod F~ 



7 ■ 



This means (T. iy'ip) = {TiyTi'\p'), where p' is the path corresponding to *(1, 0, 0, 0). In 
view of (|3.26|) . this reduced angle variable corresponds to 



















5 
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2 










1 


1 








From ^(pO = 112211112221331113213321, one gets (rj'yT}^-'(p') by a cyclic shift. The 

result yields {T^^^'Y^^^ (p) = 213321112211112221331113 in agreement with (lOSj) . Of course, 
the choice of the representative mod F-^ is not unique. 

The procedure is parallel for (rj^'^^Y^^'^ (p) with I = 1, ... ,4. The reduced angle variable of 



112211132211321133113221 = p' . 



(1)n6^(2), 



is 



(2) 



w + lOOO/i^ 
From this, one has 

(rfV°°°(p) 

^{2)\ 1000, 



/ -1 \ 

-2 
1000/ - 3 
V 997 J 



Tiy{Tiy{p"), 

15 /rp(2)\ 10 / ///\ 



mod F-y 



4/iS^)+8/if^+*(0,2,0,0), 



1 = 1, 
1 = 2, 

15/i[^^ + 10/if^ +*(1, 0,1,0), / = 3, 



lO/iJ^^ +*(0,2,0,0), 

? + 10/if ) + 
^21/ii^^ +3/if^ +*(1, 0,1,0), l = A. 



(1)n10„(2). „^ 



(1)n21/„{2)n3 



1 



(3.34) 



Here the paths p" and p'" are those corresponding to the reduced angle variables *(0, 2, 0, 0) and 
*(1,0, 1,0), respectively. They are obtained as 



2 


2 













1 












111222132111332113311122 = p" , 
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4 
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1 


1 








112211132211122133113321 = p" 



In p.34p . the evolutions under is given by 



Ti'\p") 



111222112211333111321132, 
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(T^^y(p") = 132113321122111112223311, 
(rf = 111221113221112211321333, 

(rf = 132213311113221132111221. 

Applying the cyclic shift T^^'^ to these results, one can check that ()3.34p leads to (I3.33p . 

3.4 Bethe ansatz formula from size and number of connected components 

The results in the previous section provide a beautiful interpretation of the character formula 
derived from the Bethe ansatz at g = |19j in terms of the size and number of orbits in the 
periodic An^ SCA. Let us first recall the character formula: 

(XI + • • • + = E ^(/^) -[-"^"''4'^"''-"^"'' • • • (3.35) 

A* 

, (a) , (a) T , 

^7(;.) = (detF) n ia)' ^ (e^)- (3.36) 

In (|3.35p . the sum is taken over n-tuple of Young diagrams /i = (^f(i),...,/i(")) without any 
constraint. All the quantities appearing in (|3.36p are determined from ^ by the left diagram 
in (|2.5p and (|2.6p - (|2.8p . The fact r2(^) G Z can be easily seen by expanding detF. It is vital 
that the binomial coefficient here is the extended one: 

a\ a(a - l)---(a- fe+l) 

— (a G Z, G Z>o), 



bj h\ 

which can be negative outside the range < 6 < a. Due to such negative contributions, the 
infinite sum (I3.35P cancels out except leaving the finitely many positive contributions exactly 
when L > \fi^^^\ > • • • > l^*-""^!- For example when L = 6, n = 2, coefficients of monomials in 
the expansion (I3.35P of (xi + X2 + ^3)^ are calculated as 

60x1x1x1 : (|/i(')|,|/i(2)|) = (3,1), 

n{{3), (1)) + ^^((21), (1)) + !^((111), (1)) = 6 + 36 + 18 = 60, 
Oxlx^'xl: (|/.W|,|^(2)|) = (1,3), 

0((1), (3)) + (21)) + 0((1), (111)) = 6 + (-18) + 12 = 0. 

In this way ^^(/x) gives a decomposition of the multinomial coefficients according to fi. It is 
known [191 Lemma 3.7] that ri(/i) > 1 if /i is a configuration, namely under the condition 

The formula (|3.36p originates in the Bethe ansatz for the integrable vertex model associated 
with Uq{An'^) under the periodic boundary condition. In this context, /i specifies the string 
content. Namely, fi^"^ si gnifies that there are mf^^ strings of color a and length if^^ . Under such 
a string hypothesis, the Bethe equation becomes a linear congruence equation at q = (j3.42p . 
and counting its off-diagonal solutions yields Q(/i). In this sense, the identity ()3.35p implies 
a formal completeness of the Bethe ansatz and string hypothesis at q = 0. For more details, see 
Section 13.61 especially Theorem 13.31 A parallel story is known also at q = 1 [15] as mentioned 
under Theorem 12.21 

Back to our An^ SCA, it is nothing but the integrable vertex model at q = 0, where Uq(^An^^ 
modules and row transfer matrices are effectively replaced by the crystals and time evolution 
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operators, respectively. In view of this, it is natural to link the Bethe ansatz formula ^{fJ,) with 
the notions introduced in the previous subsection like level set, torus, connected components 
(7~-orbits) and so forth. This will be done in this subsection, providing a conceptual explanation 
of the earlier observations 1251 [26l . Our main result is stated as 



Theorem 3.1. Assume the condition (|2.14p and Conjecture \'6.'6[ Then the Bethe ansatz formula 
^(n) (|3.36p counts the number of paths in the level set V{fJ.) as follows: 

IWI=J^(^) = E det^ n_ ia)] [a) ■ (3-37) 



7 



number ofT- orbits 

Here 7 = (7^°^) (^j^g"^ ^'^^ ■^wm extends over all the orders of symmetry, i.e., each 7!^^^ runs 
over all the common divisors of rnf^ and pf^ . 

The result (j3.37p uncovers the SCA meaning of the Bethe ansatz formula ()3.36p . It consists of 
the contributions from sectors specified by the order of symmetry 7. Each sector is an assembly 
of identical tori (connected components), therefore its contribution is factorized into its size and 
number (multiplicity). 

For the proof we prepare a few facts. Note that s^"^ p.2ip sending A to A' also defines an 
action of A on part alone. 



Lemma 3.1. The number of T-orbits in X^jA is i.e., the number of A- orbits in X^. 

Proof. Let y = A - x and y' = A - x' he any two elements of X-^j A. They belong to a common 
T-orbit if and only if there exists t £ T such that t ■ y = y' . It is equivalent to saying that there 
exist t £ T and a G A such that t • {a ■ x) = x' . Write x, x' as x = {uj, A) and x' = (w'. A'). 

Suppose A and A' belong to a common .A-orbit in X^. Then there exits a £ A such that 
a • A = A'. Hence we have a • x = (w, A') with some a) G X^. Since T acts on the X^ part 
transitively and leaves the X^ part untouched, there exits t £ T such that t ■ {a ■ x) = x' . 

Suppose A and A' belong to different ^-orbits. Then a ■ x and x' have different X^ parts for 
any a £ A. Hence for any t G T and a G ^, we have t ■ {a ■ x) ^ x' . 

Thus y and y' belong to a common T-orbit if and only if A and A' belong to a common 
^-orbit. The proof is completed. ■ 

Lemma 3.2. The cardinality of the set X^/A is given by 



i^^Mi= n 



|A (a)(mj )\ 



Proof. It is sufficient to show that each factor in the right hand side gives the number of 
^-orbits for each (a, i) block. So we omit all the indices below and regard ^ as a free abelian 
group generated by a single element s. 

Let A = {Xa)aez be any element of Ay{m,p). Then due to (|3.2ip we have • A = {Xa+n — 
Ai+n)oez- By (j3.23p and (|3.17p this implies that s™-/'*' . A = A and there is no < A; < m/7 such 
that s'^ • A = A. Hence the number of ^-orbits in A^[m,p) is | A^(m,p)|/(77T,/7). ■ 

Lemma 3.3. For any path p € V-y{fi), the size of each T-orbit S(p) is |S(p)| = detT^. 

Proof. This is due to Conjecture 13.31 and detTy > under the assumption (|2.14p . See the 
remark under (|3.25p . ■ 
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Now we are ready to give 



Proof of Theorem 13.11 



|'P(;u)| kJ^^—SI = y^(size of T-orbits) x (number of T-orbits) 



7 7 
Lemmas l3 . 1 1[373] ^ 



7 



detF V- rr |A / (a) Ja).| 



n(aj)G//"^i 7 (ai)eS^ 

-^^^ n (eiv(™!°'..i"')i) 

{ai)£H ^ « 

One can give an explicit formula for |A^(m,p)| by solving (I3.24p by the Mobius inversion: 

|A,Kp)|=X:^(^)(y_"/), (3.38) 

where /3 runs over all the common divisors of m and p that is a multiple of 7. Here fi is the 
Mobius function in number theory [35] defined by = = if A; is divisible by the 

square of an integer greater than one, and ^i{k) = {—ly if k is the product of j distinct primes. 
(This ^ should not be confused with the n-tuple of Young diagrams.) 

Example 3.6. Let us consider the decomposition of the level set P(/i) with /i = ((33222), (41)), 
which is the same as Examples 13.11 and 13.21 From the tables therein and det F = 4656, the 
formula (j3.36p reads 

1 /pj^^ + m[^^ - 1\ 1 /p^^^ + m^^^ - 1 



1 /p\ ' +m\' -1\ 1 / 
^7(;.) = detF^f^ 



On the other hand, ()3.38p gives 



, i\ / 4+2 i\ /r,\ /4+2 1 



1 ^ / v^/ \ 2 

4+2 



|A,(,«,4^))|=^Qj (T-"l')=^ 
|Ai(p?\m«)|=^(i) (T_"/l=36, 
|A,(pf),mf))| = |A,(pr),mV)| = 
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The formula (j3.37p consists of the two terms corresponding to the order of symmetry 7 = 
(7i^\ 72^^ 72^^) = (l) 1) 1; 1) s-iid (2, 1, 1, 1). Let us write them as 7^^ and 72, respectively. 
F~^^ = F \s given by (j3.1ip and F-y^ by (|3.27p . Their determinants are detF^^ = 4656 and 
det = 2328. Thus dOT]) reads 

\V{^j)\ = nip) = 139680 = (detF^J • 24+ (det J • 12. 

Accordingly the level set V{fj.) is decomposed into 36 tori as 

V{fi) = 24:{Z^/F^^Z^) U 12{Z^/F^^Z'^). 

The paths (13. 8p and ()3.13p in Examples 13.11 and 13.21 belong to one of the tori TL^jF-^^^ and 
Tj^/F^^U^, respectively. 



3.5 Dynamical period 



Given a path p G 'P-y(/i), the smallest positive integer satisfying (T^^ ) (p) = p is called the 

dynamical period of p under the time evolution . Here we derive an explicit formula of the 
dynamical period as a simple corollary of Conjecture 13.31 It takes the precise account of the 
symmetry specified by 7 and refines the earlier conjectures in |25| [26] for n > 2 which was 
obtained from the Bethe eigenvalues at g = 0. 

For nonzero rational numbers ri, . . . , r^, we define their least common multiple by 

LCM(ri, . . . , r,) = min(|Z n riZ n • • • n r,Z| \ {0}) . 

Given (r, /) with 1 < r < n and / > 1, define F~^\bj] (resp. F\bj]) to be the 5x5 matrix obtained 
from F^ ([3:25]1 (resp. F by replacing its (6j)th column by h^p ([5:20]) . Set 

where the LCM should be taken over only those (fej) satisfying detF[6j] 7^ 0. 

Theorem 3.2. Under Conjecture 13. 3^ the dynamical period under T^^^ is equal to M^"^^ ()3.39p 
for all the paths in V-y{^). 

Proof. Under Conjecture 13.31 the relation (rj-^^^^ (p) = p is equivalent to 

Mh\''^ = mod F^ U> (3.40) 

for any path p G V-yifi). In other words there exists z ^ ZP such that Nh^[^ = F-y z. By 

Cramer's formula, the solution z = {^f^)[bj)&H °^ equation is given by Zj'^ = J^^^^^^^^- 

The quantity A/"/*^^ (I3.39P is the smallest positive integer M that matches the condition z G 
hence ■ 

Remark 3.4. To the dynamical period under a composite time evolution T = Yiri (-^Z^^)'^' ' 
the same formula (j3.39p applies by replacing the role of h^^^ therein with ^l''^ ■ 
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Remark 3.5. The formula ()3.39p can be simplified when n = 1 (hence r = 1) and the order 
of symmetry is trivial, i.e., 7 = (7j-^^) with ^J^^^ = 1 as explained in |24^ equation (4.24)]. In 



particular if Z = 00 furthermore, one gets {p\ 
= LCM 



(1) 



0) 



l'j-l)Pgi J l<j<g^ 



which reproduces the result obtained in [lOl 



r(r) 

and compute its dynamical period A/^ under 

2328. The relevant 

vectors /ij"^^ are listed in Example 13.51 Consider T^^^ for instance. Then, h^^ = *(3, 2,0, 0) and 



Example 3.7. Let us take the path p 

the time evolution Tj^^^ . The matrix F-y is given in (|3.27p and detF^ 

i^'^^ are listed in Example 13.51 Consider Tg^^ for instance. Then, /ig 



^7 [21] 





12 


-3 


-1\ 




( 


8 


3 


-3 


-1\ 


2 


19 


-2 


-1 


, F^[12] = 




4 


2 


-2 


-1 





-6 


10 


2 




-3 





10 


2 


Vo 


-3 


2 


3^ 




V 


-1 





2 


3/ 




12 


3 


-A 




/ 


8 


12 


-3 


3\ 


4 


19 


2 


-1 






4 


19 


-2 


2 


-3 


-6 





2 




-3 


-6 


10 





V-1 


-3 





3;^ 




I 


-1 


-3 


2 


0/ 



with detF^[ll] = 840, detF^[12] = 108, detF^[22] = 276 and detF^[22] = 204. Thus (fM]) is 



calculated as 



(1) 



LCM 

LCM 



detK, 



detF-, 



detF, 



7 



detF^ 



det [11] ' det F^ [12] ' det F^ [21] ' det F^ [22] 
/2328 2328 2328 2328 



V 840 ' 108 ' 276 ' 204 



, , 97 194 194 194\ 

LCM — , , , 

'35 9 23 17 y 



194. 



Similarly, we get 



(2 



(2 



LCM 



LCM (12, 24, 24, 24) = 24, 

LCM(6, 12, 12, 12) = 12, 

582 1164 1164 1164 
17 ' 65 ' 377 
388 776 776 776 
^'Ts"' 141 '197 
291 582 582 582 
^'1T'T79'T07 
/1164 2328 2328 2328 



1164, 



LCM 
LCM 
LCM 



776, 



582, 



169 ' 49 ' 1009 ' 265 



2328. 



Despite the nontrivial order of symmetry 7 = (2, 1, 1, 1), these final results coincide with the 
trivial case 7 = (1, 1, 1, 1) treated in Example ll.il We have checked that these values agree with 
the actual dynamical periods of p by computer. 
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3.6 Relation to Bethe ansatz at qr = 

Let us quickly recall the relevant results from the Bethe ansatz at g = 0. For the precise 
definitions and statements, we refer to [19]. Consider the integrable ?7g(^n^) vertex model on 
a periodic chain of length L. If the quantum space is the L-fold tensor product of the vector 
representation, the Bethe equation takes the form: 



sm vr u. 



sm TT u. 



. n J.) r-r.^ ; (3-41) 



[uf' - \/^Mai) / j=i sin7r(n,^'''^ ~ ^ y/-ihCab) 



for 1 < a < n and 1 < i < j/x^")]. Here L > l/i*-^^] > ••• > l/U^"-*! specifies a sector (I2.20p 
preserved by row transfer matrices. The parameter h is related to g by = e"^'^''. 

Fix an n-tuple of Young diagrams, i.e., the string content fi = {fJ^^^\ ■ ■ ■ , /i^"^) as in (j2.5p . We 
keep the notations (j2.6p - (|2.9p . By string solutions we mean the ones in which the unknowns 
\ 1 < a < n,l < i < |/u("^|} are arranged as 

U U + + 1 - + ^1 \i<k<i], 

where e^^^, stands for a small deviation, uf^ is the string center of the ath string of color a and 
length i. For a generic string solution, the Bethe equation is linearized at g = into a logarithmic 
form called the string center equation: 

Aaia,bjpu^l^\{pt^ +m^^Kl) modZ (3.42) 

{bjl3)(iH 

for {aia) G H. Here the G x G coefficient matrix A = {Aaia,bj /s) {aia),{bj (5)eH is specified as 

Aaia,bjl3 = SabSij5a(3{p\"'^ + rnf"^) + Cab min {I'tKlf^) - SabSij. (3.43) 

It is known [19] that A is positive definite if is a configuration, namely if p.46p is satisfied. 

There are a number of conditions which the solutions of the string center equation ()3.42p 
are to satisfy or to be identified thereunder. First, the Bethe vector depends on u\'^^ only via 
^2iTy^ul ^ ^ Therefore the string center should be understood as u^'^ G M/Z rather than M. 

Second, the original Bethe equation ()3.4ip is symmetric with respect to n-"^ for i = 1, . . . , j/U^'^) | , 
but their permutation does not lead to a new Bethe vector. Consequently, we should regard the 
string centers of the (a, i) block as 



(a). 



Last, we prohibit the collision of string centers uj*^ = uf^^ for 1 < a 7^ /3 < m^""^ for any 

(a, i) S H. This is a remnant of the well-known constraint on the Bethe roots so that the 
associated Bethe vector does not vanish. To summarize, we consider off-diagonal solutions 
u = i'^i'^l) (^nia)&H *° string center equation (|3.42p that live in 



u; n u 



(,)) G ((M/Z)"' - A („))/© w for each (a,i) G H, 



where = {(f 1, . . . , Vm) G (M/Z)™" \ Va = vp for some 1 < a 7^ /3 < m}. For simplicity we will 
often say {q = 0) Bethe roots to mean the off-diagonal solutions to the string center equation. 
Let lA{ii) be the set of the Bethe roots having the string content /x. The following result is 
derived by counting the Bethe roots of (j3.42p by the Mobius inversion formula. 
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Theorem 3.3 ( [19\ Theorem 3.2]). Under the condition p.46p . the formula ()3.36p gives the 
number of Bethe roots, namely ^l{fJ,) = |^(^)| is valid. 

Given an extended rigged configuration r = ('^1 a) (ai)eI7 qgz ^ -^(z^)' define a map 

V t,aJ {ai)&H,aeZ ^ \"'i,a J {aia)eH 

by 



J2 ^-.,6,73^5') =^(pf)+mf) + l)+rg+a-l for (am) G ^- (3-44) 
{bjl3)€H 

Theorem 3.4. The map ^ induces the bijection between the set of angle variables and 

the setU{^) of Bethe roots. 

The proof is parallel with |24l Section 4.2] for the n = 1 case. The induced bijection will 
also be denoted by It also induces the action of the time evolutions T on ZY(/i) by Ti^\u) = 
'^(Ti^\r)^ for u = ^(r). Under the condition (12.140 . Conjecture 13.21 and Theorem 13.41 lead to 
the following commutative diagram among the level set, the angle variables and the Bethe roots: 



Vifi) Jifi) U{fi) 



T 



T 



T 



(3.45) 



P(/x) J{li) U{fi) 

Their cardinality is given by 

\v{fi)\ = \j{fi)\ = mfi)\=n{fi) 

with defined in (j3.36p . We will argue the time evolutions of the Bethe roots further in 

Section 1121 



3.7 General case 

Let V{fi) be a level set. From fi, specify the data like m\"'\ I'f^ and the vacancy number pf^'^ by 
(j2.5P " (j2.9p . In Sections l3.2H3.5l we have considered the case Vpf ^ > 1 (l2Til) . In this subsect ion 
we treat the general configuration, namely we assume 

/i = (/z(^\ . . . , /i^")) satisfies pj"^ > for all (a, i) G H. (3.46) 

It turns out that the linearization scheme remains the same provided one discards some time 
evolutions and restricts the dynamics to a subgroup T' of T. We begin by preparing some 
notations about partitions. 

Let A, j> be partitions or equivalently Young diagrams. We define A U to be the partition 
whose parts are those of A and z/, arranged in decreasing order [27]. For example, if A = (4221) 
and u = (331), then A U = (4332211). 

We denote by O^f"^ the (a, i) block in the diagram n^""^ in the sense of ([23]), and by 2n[''^ the 

corresponding block in fi^"'^ U which is a (2m^'^^) x rectangle. 

Let A be a partition. We say that A covers the block □ in ^'-"^ if D^"^^ is inside the diagram A 
when the two diagrams A and /i^") are so placed as their top-left corners coincide. 
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Definition 3.1. The block nj"^ is null if p^""^ = 0. 



(a) 



Definition 3.2. The block ^ is convex if (^("-i)u/i("+^)) covers the block 2a[''^ in /i(")u^ 
Here we interpret //'-''^ = (1^) and ^u^"^^^ = 0. 

Say that the time evolution T^^"^ is inadmissible to ^ if there exists 1 <i < Qa such that □j'^'* 
is null and convex, and if^ > I > l^^}-^. Otherwise we say that T^^"^ is admissible to /i. Here 



^i'^Li is to be understood as 0. Note that when — l)"^-, = 1, the block does not cause an 
inadmissible time evolution even if it is null and convex. We will also say that T^^'^'' is admissible 

or inadmissible to a path p S depending on whether T^"''^ is admissible or inadmissible to /i. 
As a generalization of Conjecture 12.11 we propose 

Conjecture 3.4. T^^^{J^{^)) = V{fJ,) holds if and only if T^^^ is admissible to fi. Thus, 
T^'~\v{ii)) = 'P{p) is valid for any r,l if and only if there is no block D^"^ that is null, convex 
and of length such that — > 2. 

Example 3.8. Consider the following configurations corresponding to the evolvable paths 
pj = 121122333111112222233333 (upper one) and pu = 111222133321111222233333 (lower 

(2) . (2) 

one). Hatched blocks are null and convex. They show that is inadmissible to pi, and 
and are inadmissible to pjj. 



(a) ,{a) 



(1) 



(2) 



(1) 



(2) 



Let T' be the abelian group generated by all the time evolutions xj^""^ admissible to fi. Then 
Conjecture 13.41 implies that T' acts on the level set V{fi). We define the connected component 
S'(p) to be the T'-orbit in V{i-i) that contains p. 

As it turns out, the restriction of the dynamics from T to T' will be matched by introducing 
a sub-lattice L of Z^' as follows: 



M) 



if is null and convex. 



^i+i + 1 otherwise. 



See (|3.20p for the definition of h^^\ The subgroup T' C T and here are analogues of T 
and rj^'^^ in Remark 13.31 in the present setting. It is easy to check 



(a) 
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Lemma 3.4. Let fi 



(1) 



(a, i) block of n is null and convex. Then: 



, /i^")) be a configuration, i.e., p.46p is satisfied. Suppose the 



0. 



[a) For 6 = a lb 1 there is no 1 < j < gi, such that if" > I - > l °^i- 
Recall that the matrix is defined by ()2.8p and ()3.25p . 
Proposition 3.2. is a sub-lattice ofh. 

Proof. Let Fa be the ga x {ga-i + 5a + Qa+i) sub-matrix of F-y for color a row and colors a, a it 1 

F^,hJ-i?'^'\. ... . Let / 



column indices. Its (6j)th column (6 = a, a± 1) is the array {Pai.hj hf^) i^^^g ■ '"^(a) 



. Let h^'^l E 



be the sub-vector obtained from by extracting the color a rows. We are to show that for 

any a, all the column vectors of Fa can be expressed as a linear combination of h^'^l^ with integer 

coefficients. For simplicity we demonstrate it by an example. The general case is similar. 

Suppose ga = 7 and the (a,i) block is null and convex if and only if i = 1,2,6,7. Then the 
vectors h^'^d) (1 < ^ < 7) are given by the columns of the matrix in the left hand side of the 
following equation: 





(h 


h 


k + l 


^5 + 1 k + l 


k h\ 




fh h 


1 \ 






h 


h 


k + l 


k + l k + l 


k h 




h h 


1 






h 


h 


k + l 


k + l k + l 


k h 




h h 


1 






k 


h 


k 


k + l k + l 


k h 






1 






k 


k 


h 


k k + l 


k h 






1 






k 


k 


k 


k k 


k h 






k h 






Kh 


h 


h 


h h 


h hJ 




\ 


h h) 




where li = 


'■i 


. The 


right hand side 


is obtained from the left (and vice versa 



(3.47) 



transformations adding integer multiple of one column to other columns. 

By (i) of Lemma 13.41 we have p^^^ = 0, and by {ii) of the same lemma any column of Fa is an 
integer multiple (and possible integer addition to the 4-th and 5-th rows) of one of the following 
vectors: 



h 
k 
k 
k 
k 



for I'^P 



h 
h 
k 
k 
k 



for 



(b) 



/ 
I 

X 

y 
k 



for k > I'j 



I > k 



(h\ 

h 
h 
h 
h 
h 



for 



(b) 



where x and y are some integers. Clearly these vectors can be expressed by linear combinations 
of the columns of the matrices in ()3.47p with integer coefficients. ■ 



Proposition 3.3. The volume of the unit cell of the lattice L is equal to W {l}^'^ — l^i^i), where 



the product is taken over all the null and convex blocks □ 



(a) 
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Proof. By using the above example we have 

\ (h - h 



(h 


h 


1 


h 


h 


1 


h 


h 


1 



1% Ij 

h hi 



h 1 

1 



— h h 
hi 



The general case is similar. 



Let the abelian group T' act on the torus L/F-y by (j3.30p . This action is transitive, since T' 
contains the free abelian subgroup generated by all T^|^2) 's- 

Conjecture 3.5. Under the condition ()3.46p . n V+{fJ,) ^ holds for any p G V{^). 



Owing to this property, one can introduce by the same scheme as (|3.6p with T replaced 
by T'. Similarly, let x' : T' ■ {A - x) — > L/F^Z^ be defined as in (I3.3ip . As the generalization 
of Conjecture 13.31 we have 



Conjecture 3.6. For any path p G V-y{n), the map ^'y. := x' ° ^' gives a bijection between the 
connected component T,'{p) and the torus IL/FjI^^ making the following diagram commutative: 



S'(p) 
T' 



T 



By Proposition l3.3| this conjecture implies the relation = det F^/ Y[ {li —li^i) , where 

the product is taken over all null and convex blocks nj'^^- Note that this factor is common to 
all the connected components in the level set. Thus Conjecture 13.61 also tells that 



which refines \26\ Conjecture 4.2]. 



Example 3.9. Consider the path pm = 121212343434121212 of length L = 18. Its soliton 
content is given in the following table. 



(ai) 


,W 


(a) 

m- 


(a) 
Pi 


i,a 


(a) 

Ti 






(11)* 


3 


1 








1 


3 


(3,1,0,0) 


(12)* 


1 


9 





0,... 


,0 9 


1 


(1,1,0,0) 


(21) 


3 


2 


3 


0,0 


1 


1 


(0,0,1,0) 


(31)* 


3 


1 








1 


3 


(0,0,0,1) 
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Its (1, 1), (1,2), (3, 1) blocks marked by * are null and convex. The soliton content is depicted 
as follows. (Null and convex blocks are hatched.) 

^(1) /i(') M^') 



T2^\ T^^'^ and Tg'^^ are inadmissible to pni^ while all the other T^^^^'s are admissible. The 
matrix for pjjj is given by 



/ 


6 


2 


-6 






2 


2 


-2 







-3 


-1 


15 


-3 


V 








-6 


6 J 



The cardinality of the connected component in the level set is calculated as 

\^'ipiii)\ = det FV(4^) - /«) (/« - /«) (/f ) - ) = 432/(2 • 1 • 3) = 72. 

Example 3.10. Consider the path piy = 1122331142233444 of length L = 16. Its soliton 
content is given in the following table. 



(ai) 


H 


(a) 

m- 


(a) 
Pi 




(a) 






(11) 


4 


1 








1 


3 


(3,2,0,0,0) 


(12) 


2 


4 


2 


1,1,0,0 


2 


1 


(1,1,0,0,0) 


(21)* 


4 


1 








1 


4 


(0,0,4,2,0) 


(22)* 


2 


2 





0,0 


2 


2 


(0,0,2,2,0) 


(31) 


4 


1 








1 


1 


(0,0,0,0,1) 



Its (2, 1), (2,2) blocks marked by * are null and convex. The soliton content is depicted in the 
followin 
to piv- 



(2) (2) 

following with the null and convex blocks hatched, showing that and Tg are inadmissible 



(1) 



(2) 



(3) 
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The matrix F-y for p/y is given by 



V 



-4 
-2 
8 
4 
-4 



-2 
-2 
4 
4 
-2 



\ 


-4 

-2 

8 / 



The cardinahty of the connected component is calculated as 
\Y.'{piv)\ = detF^/(/f^ - if^) [if^ - ) = 2048/(2 • 2) 

3.8 Summary of conjectures 



512. 



Here is a summary of our conjectures. Each one is based on those in the preceding lines. Thus 
the principal ones are the linearizations. 



Condition on soliton/string content /i 



Vpf ^ > 1 (imi) 



Stability of dynamics 
Intersection with highest paths 
Linearization 



Conjecture 12.11 T 
Conjecture 13.11 
Conjecture [32]~ [331 



Conjecture 13.41 T' 
Conjecture 13.51 



Conjecture 



Relevant torus 



4 Tropical Riemann theta from combinatorial Bethe ansatz 

4.1 From tropical tau function to tropical Riemann theta function 

Let (^'-"^ r^")) be the color a part of the rigged configuration (/x, r) depicted in the right diagram 
of (12. 5p . We keep the notations in p.6p - p.9p . An explicit formula for the image path of the 
KKR map (p~^ (j2.18p is known in terms of the tropical tau function [23]. It is related to the 
charge function on rigged configurations and is actually obtained from the tau function in the 
KP hierarchy [11] by the ultradiscretization with an elaborate adjustment of parameters from 
the KP and rigged configurations: 



Tk,d 



minji Y: C;,min(t),/f)ivgiv(2+ J^rg 



) ^.(«) 

a ta 



(aia), {bjl3) 



(aia) 

[l<d<n + l) 



(4.1) 



(ia) {jP) ) 

Here the sums range over (aia) G H wherever N^'^^ is involved. Thus in the second line of (|4.ip . 



(ia) runs over l<i<(?i,l<a< m[^^ and so does (j/3) over 1 < j < g^, 1 < /3 < m^"''. The 

last term X^q-^) l^f^ N^'^^ is to be understood as when d = n + 1. In ()4.ip . min is taken over 

-^ici ^ {0' -*-} fo'^ {aia) G H. Thus it consists of 2^ candidates. Obviously, Tk,d £ ^>o holds. 

In the present case, the paths are taken from [B^'^)®^ . We parameterize the set B^'^ = 
{1,2, . . . ,n + 1} as {(xi,X2, . . .,Xn+i) G {0, 1}"+^ \xl^ h = 1}. 

Theorem 4.1 (|23|). The image of the KKR map cj)^^ is expressed as follows: 
(/>"^ : (/i, r) xi X2 • • • XL e (S^'^)®^, 

Xk = {xk,l,Xk,2, ■ ■ ■ , Xk,n+l) G B^'^, 
Xk,a — Tk,a Tk—l,a Tfe a— 1 ~l~ 'Tk—l,a—l 



(d) 



(2 < a < n + 1). 



(4.2) 
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In the context of the box-ball system, Xk,a represents the number of balls of color a in the 
A;th box from the left for 2 < a < n + 1. The remaining x^^i is determined from this by 
Xk,i = 1 — Xk^2 — ■ ■ ■ — Xk,n+i, which also takes values x^^i = 0, 1 in the present case. 

The following is a special case (VA^ = 1) of [23l Proposition 5.1]. 

Proposition 4.1 (|23]). The tropical tau function satisfies the tropical Hirota equation: 

Tk-i,d + Tk,d-i = max(rfc,d + Tk-i,d-i, Tk,d + Tk-i,d-i - 1) (2 < d < n + 1), 
where Tk,d is obtained from Tk^d by replacing r^^^ with r^^^ + /-"^(Jai- 

Lemma 4.1. Let p,p' € (^B^'^^^^ be the highest paths whose rigged configurations are r) 
and (/i',r'), respectively. Then the rigged configuration of the highest path p iSSi p' G ^B^'^^^"^^ is 
(/i,r) U {fi',r'), where (f') = (r"!*^' Pi"^^) (aia)eH '^'^'^ pf^ vacancy number of {ii,r). 

In the lemma, U stands for the sum (union) as the multisets of strings, namely, the rows 
of ^ attached with riggings. More formally a string is a triple {o.jl^i^\fYa) consisting of its 
color a, length l^f^ and rigging r'f^ that is labeled with H. Let p G (S^'^)®^ be the highest path 
corresponding to the rigged configuration {n,r). From Lemma |4.H the rigged configuration of 
the highest path p®^^ G (^^j.^) U (/i, r^) u • • • U (/i,r*^), where r'= = (rj^j + (/c - 

^)Pi"^^) (aia)&H' Pictorially, this corresponds to extending the right diagram in (12. 5p vertically M 

times adding pf^\2pf^\ . . . , (M — l)p["'^ to the riggings in the (a, i) block for each {a,i). This 
reminds us of (|3.ip . and is in fact the origin of the extended rigged configuration in Section [3.11 
We proceed to the calculation of the tropical tau function r^'^ associated with the above 

rigged configuration (;U, r-^) U r^) U • • • U r^^) . In (j4.ip . the variable A''^'-^'' is to be replaced 

by the M replicas N^'^^^, . . . , N'f'^ to cope with the M-fold extension: 

r^,d = —^A\ E E C^^{l?.lf)NQ^,Nfi^ 

I (aid), {bjli) l<s,t<M 



a,s 



(aia) l<s<M 



^E E <. + E E f (4-3) 

Ha) l<s<M as) l<t<M ) 



(ia) l<s<M 

All the summands here are invariant under permutations within N-\^-^, . . . , N-^^ for each (a, i) 
except the replica symmetry breaking term "^(^^ia) J2i<s<Mi^ 

- l)pf^N^'ts- Due to pf^ > 0, 

^le{o,i}to"ng 



Jf^N^;^^, Due to ) 

one can reduce the minimizing variables N^"^ ^ G {0, 1} to nl"^ G {0, 1, . . . , M - 1} such that 



^(a) ^(a) ^ . . . ^ jy{a) ^ ^(a) = jyi'') =n (44) 



As the result, ()4.3p becomes 

^A/ _ I 1 ^ /r(a) Ah)\ (a) [b] 

^fe,d--™™<^- 2^ Cab mm (^/- j^-i.Q^j,/? 



(a) / (a) 

{aia) \ / (la) 



, Sr^ ( (a) (a) . ^tli'^tl - {a)\ , (1) , (rf) 
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where the minimum is now taken over nf^^ G {0, 1, ... ,M — 1} for all the blocks {a,i). The 

notation can be eased considerably by introducing a quadratic form of n = {v!f^ 
follows: 



{aia)&H 



as 



B 



min <J — *ni?n + ( r 



2 '^"-1 ^ "oo 



{aia),{bj(S)&H^ 



P 



{l<d<n + l,l<k< ML), (4.5) 

(4.6) 
(4.7) 



h 



'5-min(/,/f)))(^^^)^^ (l<c<n,^>l), h 



(n+l) 



0. 



(4. 



(c) 

The G-dimensional vector here should be distinguished from the (^-dimensional one 



h^^^ (j3.20p . It is the velocity vector of the time evolution T^'^' in i7(/u) and i7(/i) in the light 
of (|3.3p . The matrix B is symmetric and positive definite. See for example [T9t Lemma 3.8]. 
We call B (|4.6p the tropical period matrix although a connection to a tropical curve is yet to be 
clarified. We also introduce the G-dimensional vector 1 = (l)(aio)e-ff • Then from (|2.7p and ()4.8p 
we get 



Bl = 

-^'■'■rv^ 



( E (5a.5.,5./3pf^ + a,min(/f),/f))') 



aia) 



Lh 



(1) 
1 ' 



(aia) 



(4.9) 



(4.10) 



Example 4.1. With the ordering of indices in H (j3.10p . the tropical period matrix B (|4.6p and 



the vectors 


r, P 


(021) 


for the rig 


ged confi^ 


juration (13. 9p are given as follows: 




/lO 


6 


4 


4 


4 


-3 


-1\ 




/4\ 




/4\ 




6 


10 


4 


4 


4 


-3 


-1 




2 




4 




4 


4 


11 


4 


4 


-2 


-1 




6 




7 


5 = 


4 


4 


4 


11 


4 


-2 


-1 


r = 


5 


P = 


7 




4 


4 


4 


4 


11 


-2 


-1 




1 




7 




-3 


-3 


-2 


-2 


-2 


10 


2 









2 




V-1 


-1 


-1 


-1 


-1 


2 


3^ 




\o) 







(4.11) 



The velocity vector h^""^ ()4.8p for the time evolution Tl^""' is specified as 



(a) 



hi 



(1) 



(1) 



.(2) 



-h 







2 



/1\ 

1 
1 
1 

1 





n>3 



/3\ 

3 

2 
2 
2 


Vo/ 



(2) 



.(2) 







3 

VJ 



hi 








4 

vv 








1 
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Now we take M to be even and make the shifts n — t- n + ^1 and k ^ k' = k + 
Using (|4.9p . we find 

Tfcf = - min I ^'nBn + * - | - fchS'^ + h^;^)) n| + um 'h^^h + VMk' + wm- (4.12) 

The first term is formally identical with (j4.5p but now the minimum extends over —^<nf^^<^. 
The scalars um, vm, wm are independent of k' and d, therefore these terms are irrelevant when 
taking the double difference r^^^ - r^Li^a ~ "^k'^a-i + '^k'^-i,a-i ^ ™ <\4:.2\\ . 
Passing to the limit M — >• cxd, we see that the relevant term in (I4.12p tends to 

hm (t,¥, - UM *hWl _ y^k' -yjM)=e(r-^- kh^l^ + hW) . (4.13) 

M: even— >oo \ ' / \ 2 / 

Here denotes the tropical Riemann theta function 

e(z) = - min I -*nBn + *znl , (4.14) 
nez<3 [2 J 

which enjoys the quasi-periodicity 

e(z + v) = e(z) +*v5-^ (^z + I) forvG5Z^. (4.15) 

In the context of the periodic box-ball system, the tropical Riemann theta function was firstly 
obtained in this way in |20] for rank n = 1 case. See (9l [TOl [30] for an account from the tropical 
geometry point of view. Another remark is that Proposition 14.11 and (|4.13p directly lead to the 
tropical Hirota equation for our 0: 

e(j + hW + h(^-i)) + e(j + hi^) + hW) 

= max{e(j + h« + hW) + e(j + + h(^"i)), 
e(j + hS^) + h« + h(^-i)) + e(j + hW) _ 1}, 

where we have set J = r — ^ — /ch^^\ For n = 1, see also [10| . 
From (14.21). we arrive at our main formula in this section. 



Theorem 4.2. The highest path xi (8) X2 'S' • • • (8> xl G V+{fi) corresponding to the rigged con- 
figuration {fi, r) is given by Xk = (xk^i, Xk,2, ■■■ , Xk,n+i) with Xk,i H h Xk,n+i = 1 and 

Xk,a = (r - f - kh^i^ + h^S^) - e (r - I - (fc - l)h« + h^^^) (4.16) 
_ e (^r _ I _ A;hS') + h(2"i)) + e (r - I - (fc - l)hS') + h(^-i)) (2 < a < n + 1). 

The data r, p and h^*^^ are determined from (/i,r) by ()2.5p . (j2.7p . ()4.7p and (j4.8p . 

Compared with the previous expression (|4.2p with (j4.ip . the alternative formula (|4.16p en- 
joys an extra symmetry under the shifts r — t- r + BI^'^ and k ^ k + LZ due to the quasi- 
periodicity and ([O]) . In other words, (|i36]) is defined for r G iP / BZ^ and A: G Z/LZ. 
In view of this, it is natural to relate Theorem 14.21 with our periodic A^^ SCA. 

Theorem 4.3. Assume the condition (j2.14p . Letp = xi®X2®- ■ - ^xl G V{^) he the path of the 
periodic A^'' SCA whose action variable is the quasi-periodic extension (|3.ip ofr= {'''{"(i) (aia)&H ' 
Under Coniecture \2>.2\ Xk = {xk,i,Xk^2^ ■ ■ ■ ,Xk^n+i) is given by (j4.16p . 
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Note that this reduces a solution of the initial value problem p — )• Tj^^^'^ ■ ■ ■ T^^J^^ (p) to a simple 
substitution r ^ r + hf'''^ H h h\''''\ 

Remark 4.1. In addition to r — )■ r + BlP , each term in (|4.16p is invariant under r = 
(^>)(a*a)eH ^ (^>(a))(aia)e// ^^^^ permutation a G 6^(.) that can depend on {ai) G 
Here <Sm denotes the symmetric group of degree m. In fact the set i7(/x) of angle variables is 
naturally described as 

J{p.)^lPlBlPl n_®„^W- (4-17) 

(ai)Gl7 

To go to the right hand side, one just forgets the inequality r^*^ < v'f^j^^ within each block 



{ai) G H identifying all the re-orderings. Compare also the matrix B (j4.6p with ([37 

When n = 1, one can remove the assumption 'Under Conjecture 13. 2f in Theorem 14.31 due 
to [21]. If further Vm^-^^ = 1, the formula (|4.16p coincides with [20l equation (3.8)]. Theorem l4.3l 



here tells that it also holds throughout > 1 by a natural extrapolation whose concrete form 
is given in (|T26]) and (jOTl) . 

The case n = 1 and general m|^^ in ()4.16p takes a different form from the corresponding 
formula in [21] expressed by a higher characteristic tropical Riemann theta function. This 
apparent difference in guises is caused by a different choices of effective minimizing variables 
as nfl^ in (|4.4p . In fact, introducing the M replicas -^f i j • • • j -^f Af ^ ™ (|4.3p is not the 
unique way of handling the tropical tau function for the M-fold extended rigged configuration 
(^, r^) U r^) U- ••□(/!, r*^) . Another natural option is to stick to the form (14. ip but extend A''^^^'* 

from 1 < a < nif^ to 1 < a < Mmf^ assuming (13. ip for the rigging. Then the minimizing 
variable n^"^ can be introduced simply via 

At/^) = Ar(") = . . . = A^*^) , = 1, A^(") , = • • • = AT^-^) = 0. 

' ' i.n. ' %jn. '+1 I, Aim- ' 

What makes the calculation tedious after this is that one has to classify nf^ by modm^"\ 



This was done in |21j for n = \. Here we omit the detail and only mention that the result is 
expressed in terms of a rational characteristic tropical Riemann theta function with the g x g 
reduced period matrix: 

B-' = <^^g{mt\,^,Tj F = {5^,5., Jf^j^ + C., min (/f \ /f ^-f ) ia.),m,H- 

For n = 1, this coincides with |2H equation (5.1)]. 

Let us proceed to an explicit formula for a carrier of type B^'K Consider the highest path 
p = xii5<)X2'S<)---'S'XL G V+{pi) in Theorem 14.21 expressed as (|4.16p . We consider the calculation 
of the time evolution tI^\p) = p' = x'l • • • x'^ according to (12. 2p . Locally it is depicted as 



Xk Xk+l 



Vk-i 



Vk 



^'k ^'k+1 



From the proof of Lemma [D.H a carrier satisfying the periodic boundary condition yo = Ul ^ 
B^'^ and thereby inducing the time evolution T^^^ can be constructed by n^'' ®p — p®yQ. Here 
^1.' g ig lY^Q semistandard tableau of length I row shape whose entries are all 1. This fixes the 
carriers in the intermediate stage yi,y2, ■ ■ ■ , yL-i G -B^'' by yo^{xi^- ■ -^Xk) i— )• {x\^- ■ ■®x'i^)®yk 
under the isomorphism B^'^ ig) (^B^'^^®'^ ~ ^^i.i)'^'^ ^ B^'K 
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Theorem 4.4. The carrier G -B^'' corresponding to |^.i6p in the above sense is given by 
Vk = {yk,i,yk,2, yk,n+i) with y^^i -\ \- yk,n+i = I and 



= e (r - B _ kh? + h(^)) - e (r - H _ feh« + + h(2) 



(4.18) 

_ e (^r - I - fchS') + h(^-i)) + G (r - I - A:hi') + hf'^ + h^^-^^) (2 < a < n + 1). 

This can be verified by modifying the derivation of Theorem 14.21 slightlv by using \23\ Theo- 
rem 2.1] and |18| Lemma 8.5]. The periodicity yo = y^ is easily checked by ()4.15p . Natu- 
rally ()4.18p gives the carrier for general paths in Theorem 14.31 under Conjecture 13.21 For n = 1 
such a result for a higher spin case was obtained in [22] under the condition Vm^^^ = 1. We will 
present a modest application of the formula (j4.18p in Section 14.41 



4.2 Bethe vector at q 



We assume the condition (|2.14p and Conjecture 
equation (|3.44p simply as 



from tropical Riemann theta function 

in this subsection. Write the string center 



Au = c + r + p (4.19) 

using the G-dimensional vectors c = {^{pf"^ + m,-"'* + and p = {a — l)(aia)£H- Note 

that the time evolution of the angle variable ()3.3p is written as T^^\r) = r + h['^'' with h 



defined by (|4.8p . In view of (I4.19p . the time evolution of the Bethe roots introduced in (13.45P 
is expressed as Ti'^\u) = u + A~^lif\ At first sight, this appears contradictory, because T^^^^ 
is a transfer matrix at g = 0, which should leave the q = Bethe eigenvectors invariant up to 
an overall scalar hence the relevant Bethe roots as well. The answer to this puzzle is that the 
path p G Vl^) that we are associating with u or r by p = ^~^{fi, r) is a monomial in (C""*"^)®^, 
which is not a Bethe vector at g = in general. 

It is easy to remedy this. In fact, for each Bethe root u or equivalently the angle variable 
r = — c — p G '7{fJ'), one can construct a vector |r) G (C"^-^)'®^ that possesses every aspect 
as a (/ = Bethe vector as follows: 



seJip.) 



v(s). 



v(s) 



exp I —27ry 
( 3^1, i(s) 

\xi,„+i(s)^ 



T*s 



A Xr + c + p) 



(4.20) 
(4.21) 



Here 1 = (l)(aia)e/i' in (14.21|) . which has also appeared in (14. 9p . The component Xfc^a(s) G {0, 1} 
is specified by Theorem 14. 3[ Namely, 



Xk,a{^) 



e IS 
-eis 



p 



p 



2 A;h(^)+h(s))-G(s 2 



P 



{k 
P 



2 A;hS^)+h(S-i))+G(s 2 



l)h«+h(^-i^ 



for 2 < a < n + 1 and Xfc^i(s) + • • • + xj^^^+ii^) = 1- The vectors p and h['''* are defined in ()4.7p 
and (|4.8p . We embed V{^) into (C^+i)®^ naturally and extend the time evolutions to the latter 
by C-linearity. The vector v(s) here is the path corresponding to the angle variable s G J{^j) 
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in TheoremSSl It follows that r/''^(v(s)) = v(s + hj*^^). Then using J{ij) + = J{n), it is 
elementary to check 



T/^)(|r))=c_^„^Jr), (4.22) 



exp (^2TrV^ ^hf") (^A-^r + c + p) + ^ = exp (^IttV^ *hf'^ + ^ . 



(r) 

The last expression of the quantity c (r) reproduces the q = Bethe eigenvalue A, = 

-hj , r 

Yija ( ~ -ZjQ )™'°''"''''* given in [26l equation (3)] with the string center zj^^ = exp (^2ir\/—lu^jJ^) . 
Note further that the transition relation (j4.20p is inverted as 

where Cs,r denotes the complex conjugate of Cs,r- It follows that the space of the 'g = Bethe 
vectors' |r) coincides with the space spanned by the evolvable paths of the prescribed soliton 
content /i. 

C|r) = Cp. (4.23) 

Thus the approach here bypasses a formidable task of computing the g — )■ limit of the Bethe 
vectors directly, and leads to the joint eigenvectors |r) of the commuting family of time evo- 
lutions. They form a basis of the space having the prescribed soliton content and possess the 
eigenvalues anticipated from the Bethe ansatz at g = 0. These features together with Theo- 
rem 12.21 constitute the quantitative background of the identity (j2.19p . 



4.3 Bethe eigenvalue at g = and dynamical period 

Let us remark on the relation of the Bethe eigenvalue (|4.22p to the dynamical period. In [25^ 126]. 
the formula ()3.39p for the dynamical period of periodic An^ SCA with n > 1 without the order 
of symmetry 7 was found by demanding the condition 

(c_hW,r)^''" = ±l- (4-24) 

In fact, from the middle expression in (|4.22p this condition is satisfied if 

AA/"\(")=0 modAZ^. 

By the same argument as the proof of Theorem 13.21 one finds that the smallest positive integer 
satisfying this condition is 

where A[bjf3] denotes the G x G matrix obtained from A (|3.43p by replacing its {bjf3)th column 
by hf''^ (jM]). The LCM in the first (second) expression should be taken over only those {bjf3) 
such that detA[6j/3] 7^ ((6j) such that det F[bj] 7^ 0). The second equality is due to the 
structures of matrices A and F ()2.8p and verified by a direct calculation. The expression A/"/*^^' 
derived in the heuristic approach \25\ captures the main structure of the full formula (|3.39p . 
However it neither fixes the sign in (j4.24p nor takes the order of symmetry 7 into account. 
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These shortcoming are fixed of course by refining the construction of joint eigenvectors along 
the connected components. Instead of trying to split the sum (|4.2U|) into them, we simply 
introduce 

|a;,A)= Yl exp(-2^V^V^-7'('^ + dA))v((A,A) G (C"+i)®^ (4.25) 

for each angle variable (w, A) G X^/A having the order of symmetry 7. Here v(<^, A) € V'y{^) 
is the path corresponding to (</>, A) G X~^/A, and d\ G can be chosen arbitrarily. Then by 

noting T^^^w{(l), A) = w{<j) + h!f \ A), it is straightforward to check 

tI'\\lo, A) ) = exp [2^^/^%f^ F-\u + dx))\to, A), 

This eigenvalue is indeed an Mi '^\h root of unity due to the proof of Theorem 13.21 

Given A, there are det-F^ independent vectors |w, A) (j4.25p . On the other hand, the number 
of the choices of A G X^/A is given by Lemma 13.21 Obviously these vectors are all independent 
because the set of monomials involved in \uj, A) and \u}' , A') are the set of paths that are disjoint 
if A 7^ A'. From this fact and (j3.37p . we obtain the refinement of (j4.23p according to the order 
of symmetry 7 and further (bit tautologically) according to the connected components: 

C|^,A)= Cp, C|a;,A)= Cp, 

where for example po = ^~^{{uj = 0, A)) G V-y{iJ.). 

4.4 Miscellaneous calculation of time average 

As a modest application of the formulas by the tropical Riemann theta function, we first illus- 
trate a calculation of some time average along the periodic box-ball system (n = 1). Analogous 
results will be stated for general n in the end. We use the terminology in the periodic box-ball 
system. 

Let p G (i?^'i)®^ be the path with the angle variable I. The time evolution T^^'^ will simply be 
denoted by T/. SetT/(p) = xi{t)<Si- ■ -i^XLit). For n = 1, one can label Xfc(i) = {xk,i{t), Xk,2{t)) G 
B^'^ just by a;^ 2(^)5 which from now on will simply be denoted by Xk(t){= 0, 1). It represents 
the number of balls in the kth. box. Then (j4.16p reads 

xkit) = 0(j - khi + thi) - e(j -{k- i)hi + thi) 

- 0(J - khi + thi + hoo) + e(J - {k - l)hi + thi + hoc) (4.26) 

with J = I — ^. The notations (j4.6p - (|4.8p are simplified hereafter as 

B = {SijPi + 2min(Zi,/j))(jc,),(j/3)e//, H = {{i,a) \ 1 <i < g,l < a < rm}, 

9 

P = {Pi){ia)eH, Pi = L -2Ymm{li,lj)mj, h; = (min(/, /j))(j„)g^. (4.27) 

i=i 

Here li, rrii are the shorthand of rni'^ that specify the action variable (single partition) 
/i = /i^^) as in ()2.5p . The relation (j4.9p reads 



Bhi = Lhi. 



(4.28) 
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Note that the total number of balls at any time is \fi\ = Yli=i h^i- This also follows immediately 
from (|4.26!) as 

L 

xk{t) = G(j - Lhi + thi) - e(j + thi) - e(j - Lhi + thi + hoo) + e(j + thi + hoc) 

k=l 

= -L'hiB-^ (^J + thi - ^) + L'hiB^' (^J + thi + hoo - 
= L*hiB~^hoo = *hihoo = 1^1, 
where we have used (j4.28p and the quasi-periodicity (j4.15p . We introduce the density of balls: 

Let Ml be the dynamical period under Ti. Thus Xk{Mi) = Xfc(O) holds for any 1 < k < L. 
Proposition 4.2. The time average of Xk{t) under Tex;, over the period Moo is given by 

Xkit) = P- (4.29) 

■'^^ t=o 

Proof. Using (j4.26p with I = oo, we find that Yl^^~^ Xk{t) is equal to 

e(J - A:hi) - e(J - khi + AAoohoo) - e(J -{k- l)hi) + e(J - (A; - l)hi + AAoohoo). 

Note that the assumption implies Afoohoo £ BlP . Thus one can reduce this by applying the 
quasi-periodicity (j4.15p . obtaining A/'oo*hoo-B~"^hi =MooP- B 

Actually (|4.29p follows at once without this sort of calculation if the left hand side is assumed 
to be independent of k. However, such homogeneity under spatial translation is not always 
valid for some time evolution T/ with I < oo and the initial condition that possess a special 
commensurability. 

Denote the time average Yl't^o^ Qi^) °f ^ quantity Q(t) under Ti by {Q)i. Then a trivial 
corollary of Proposition 14.21 is 



pL UJ 



In particular, this average vanishes if a; is a nontrivial Lth root of unity. 

Let us proceed to a less trivial example. In [22], the number of balls in the carrier for any time 
evolution Ti was expressed also in terms of the tropical Riemann theta function. The carrier that 
induces the time evolution p i— )■ Ti{p) is the element vi E B^'^ such that vi^^p c:t Ti{p) vi under 
the isomorphism of crystals S^''® (i?^'^)*^^ ~ See the explanation around (j2.2p . 

There uniquely exists such vi for any p with density p < 1/2 as shown in \24:\ Proposition 2.1]. 
Let us consider the intermediate stage of sending the carrier yo{t) := vi to the right by repeated 
applications of combinatorial R: 

vi xi{t) (g) ■ ■ ■ Xk{t) Xk+i{t) • • • (8) xiit) 

~ + 1) (g) • • • (g) + 1) (g) yk{t) (g) Xk+i{t) • • • (g XL{t). 

This is depicted locally as 
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■■■ Vkit) ••• 

Xk{t+1) Xk + l{t + l) 

We identify yk{t) = (yfc,i(i), yfc,2(i)) S B^''' with the number of balls in the capacity / carrier. 
Namely, yk,2{t) G {0, 1, . . . , /} wih simply be denoted by yk{t). (Hence yk,i{t) = I - yk{t)-) Then 
for Xk{t) given in (|4.26p . one has [22] 

yk{t) = e(J - fchi + thi) - e(J - A;hi + (t + l)h;) 

- e(J - fchi + th.1 + hoo) + e(J - fchi + {t + l)h; + hoo). 

The preceding result (j4.18p reduces to this upon setting n = 1, a = 2 and r — ^ = J + thi. 

Proposition 4.3. The time average of the number of halls in the carrier with capacity I at 
position k under Ti is given by 

{yk)i = %B-^h^. 

Proof. The proof is parallel with the one for Proposition 14.21 by using Mihi e BZ'^ and the 
quasi-periodicity (j4.15p . ■ 

Note that the result is independent of the order of symmetry 7 as well as k. 

Example 4.2. Take the path p = 111222221111222222221111111111111112222111211 of length 
L = 45 from [21^ Example 5.2]. The action variable is 



The data from (14.270 read 

/min(/,9)\ 
min(/, 4) 
min(/, 4) 
\min(/, 1) / 

We list average {yk)i with the dynamical period A//. 



/ 


1 


2 


3 


4 


5 


6 


7 


8 


00 


M 


45 


1665 


333 


1665 


1665 


31635 


3515 


31635 


6327 


{yk)i 


p=| = 0.4 


147 
185 


44 
37 


293 
185 


6646 
3515 


1545 
703 


8804 
3515 


9883 
3515 


S = 3.11863... 



Now it requires little explanation to present the analogous result for An case. Denote 
the time average Yl't^o ^ Qi^) °1 ^ quantity Q(t) under the time evolution T^^"^ by {Q)i. 

Here A/"/^^ is the dynamical period under T^^^ . Then, the same calculation as Proposition 14.31 
using (j4.18p leads to 

Proposition 4.4. Under Conjecture 13. 2|. the time average of the number of halls of color a in 
the carrier of type i?^'' at position k is given by 

{yk,a)i = *hp)s-i (h^^-i) - h(2)) (2 < a < n + 1). 



B 



/27 8 8 2 \ 

8 27 8 2 

8 8 27 2 

\2 2 2 39/ 



P 



/9\ 

19 

19 
\37/ 
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Again this depends neither on k nor on the order of symmetry 7. It reduces to Proposition l4.3l 



for n = 1 due to h 



(n+l) 



0. See ()4.8p . In particular at Z = 1 we have 



{yk,a) 



1 



{a-l)| 



1/^ 



by means of (j4.9p and (j4.10p . From (|2.20p . this is the density of the color a balls in the 
path (i?^'^)*^^, which is natural in view of (j2.3p . The average of the total number of balls 
)l = n, D iioo , wnicn aRam reauces to jrroposition lor n = 1. 



is Yltl{yk,a)i = *h}^^S-^h^\ which again reduces to Proposition 14.31 for 



Example 4.3. Consider the paths in Example 11.11 (n = 2). The configuration (i is given in 
Example 13.21 which is the same as the one in Example 13.11 The matrix B and h["^ are listed 

in Example 14. H which says hg = . Denote the density by pa = — j^-^ — - ■ The average 

{yk,a)i is shown in the following table with the dynamical period given in the end of Example 1 1.1[ 





1 


2 


00 




24 


12 


194 


{yk,2)i 


P2 = ^ = 0.291 . . . 


^ = 0.583. . . 


194 u. (»o... 


{yk,3)i 


P3 = ^ = 0.208... 


^ = 0.416... 


li = 0.561... 



A Row and column insertions 

Let T be a semistandard tableau. The row insertion of the number r into T is denoted by T r 
and defined recursively as follows: 



T 



r = 



^l^2 



T 



if ii < r 
finished 



T 



^k 



if ik-i <r<ik 
to be continued 



Here we have denoted the entries in the first row of T by ii < 12 < • • • ^ ii, and the other part 
of T by T. The finished case includes the situation T = 0. In the to be continued case, r bumps 
out the smallest number that is larger than r (row bumping). The tableau T r is obtained 
by repeating the row bumping until finished. 

The column insertion of the number r into T is denoted by r — ?■ T and defined recursively as 
follows: 



T 



T 



Ik 



T 



r 



\i il < r 
finished 



if ik-i <r <ik 
to be continued 
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Here we have denoted the entries of the first column of T by ii < ^2 < • • • < and the other 
part of T by T. The finished case includes the situation T = 0. In the to be continued case, 
r bumps out the smallest number that is not less than r (column bumping). The tableau 
r — )• T is obtained by repeating the column bumping until finished. 



B Proof of Theorem UA\ 

In this appendix, we use the notions from crystal theory. See [14^ [131 El] ^-nd [241 Section 2.3] 
for the notations ^i,ei, the Kashiwara operator and the Weyl group simple reflection Si. The 
crystal B^'^ can be parameterized as B^'^ = {u = {ui, . . . , G {Z>q)"-~^^ \ ui+- • ■+Un+i = I}, 

where Ui denotes the number of letter i contained in a length / row shape semistandard tableau u 
[3HI8]. The cyclic automorphism a acts as o"((ni, . . . , Un+i)) = {u2, ■ ■ ■ , Un+i,ui), or equivalently 
Ui I— 7- itj+i. The automorphism a acts on a tensor product component-wise. The indices are to 
be considered in Z/(n + 1)Z. Recall that B = (S^'^)®^ and Bi is defined in ([2T5|) . We invoke 
the following factorization property of the combinatorial R into simple reflections. 

Theorem B.l ([8l Theorem 2]). For any path p £ B and u = {ui,U2, ■ ■ ■ ,Un+i) G B^'^ with 
ui ^ U2, ■ ■ ■ , Un+i, let B^'^ (g)i?9u(8)p~p(8)nGi?(E) B^'^ be the image of the isomorphism. 
Then u®p is expressed by a product of simple reflections as 

U®p= (cr® a)S2 ■ ■ ■ SnSn+l{u(g) p). (B.l) 

Of course the actual image p^u is obtained by swapping the two tensor components of (IB.ip . 
For / sufficiently large, fix the element u = [ui, . . . , Un+i) G i?^'' according to p G i?i by 

Ui = ipi{Si+iSi+2 ■ ■ ■ Sn+i{p)) (2 < i < n + 1), 

ui = l - {u2-\ h Un+l). (B.2) 

This is possible, i.e., ui > 0, if / is sufficiently large. We are going to show that substitution 
of dEH) into (IRT]1 leads to 

u = u, p = aS2S3- ■ ■ Sn+iip). (B.3) 

Since p is Tl^^-evolvable by the assumption of Theorem 12.1"! u = u means that u ()R2]l IS a proper 
carrier and 

T«(p) = a5253---5„+i(p). (B.4) 

Let Tj be the Weyl group simple reflection acting on B component- wise. Then, by using r? = id 
and crr2r3 • • • r„+i = id, the formula ()B.4p is rewritten as 

T^\P) = K2KS • • • Kn+l{p), Ka = rn+l ■ ■ ■ ra+iraSaVa+l ' ' ' r„r„+i. 

It is straightforward to check that the gauge transformed simple reflection Ka is equal to Ka 
defined by the procedure {i)-{iv) explained before Theorem 12.11 by using the signature rule 
explained in [Ml Section 2.3]. 

It remains to verify ()B.3p . We illustrate the calculation for n = 3, for the general case is 
completely parallel. Then (IB.2P reads (note S4 = 5'o,(/34 = ^0,64, = So, etc.) 

U4: = ^0{P), U3 = ^3{So{p)), U2 = ^2{S3So{p)). 

We set 



(^p(i) = 5o(tt(8)p), = 535o(n®p), u'-^^ ^ p'-^^ = S2S3So{u^ p), 
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04 = eo (p) , "3 = £3 (p^^^ ) , a2 = e2 (p^^^ ) . 

We know u\ is sufficiently large and from p £ Si, 04 > U4 is valid. To first compute Sq{u®p), 
we consider the 0-signature: 

u®p= {ui,U2,U3,Ui)®p 0: 
Thus Sq actually acts as Sq = Eq^""" ® e^^"^* = Bq^""'* Sq, leading to 

Ml U3 as U3 

^ pW = (^U4, U2,U3,ui) (g) Soip) 3 : ' • ' H~+T'^~~^^ H~^+, 

where we have depicted the 3-signature on account of 9^3 (p^^^) = ip3{So{p)) = U3. From p € Bi 
we see that 03 = £3{p^^^) = £3{Sq{p)) > </?3(5'o(p)) = ^3. Thus the next 5*3 actually acts as 
53 = e^i-"3 e^3-"^ = e^i-"^ ® ^3, leading to 



^(2) ^ p{2) 



{ui,U2,ui,U3) S3So{p) 2: h + H h + 



where we have depicted the 2-signature on account of 9J2(p^^^) = f2{S3So{p)) = U2- From p & Bi 
we see that a2 = £2{p^'^^) = £2{S3So{p)) > 'P2{S3So{p)) = U2- Thus the next 5*2 actually acts as 
52 = e^'-""^ O e^^-"^ = e^^-"^ O ^2, leading to 

^(3) 0p(3) = (n4,ni,n2,n3) S2S3So{p). 

Thus (fB?T]l gives 

U®p = (t{{u4,Ui,U2,U3)) (g) (7S'2S'35'o(p) = (ui , U2 , U3 , U4^) (g) aS2S3So{p), 

yielding (|R3]l . 



C KKR bijection 

C.l General remarks 

The original KKR bijection [T6] is the one between rigged configurations and Littlewood- 
Richardson tableaux. Its ultimate generalization in type An^ corresponding to B^"^'^'^®- ■ ■0B^^'''^ 
is available in [T^. In the simple setting of this paper, the Littlewood-Richardson tableaux are 
in one to one correspondence with highest paths in B = (i?^'^)®^. The KKR bijection in this 
paper means the one ()2.18p between rigged configurations and those highest paths. See [231 
Appendix C] for an exposition in a slightly more general setting B^'^^ ■ ■ ■ B^'^^ . 

Recall that a rigged configuration (/i, r) is a multiset of strings, where a string is a triple 
(a,/|"\r|^) consisting of color, length and rigging. See the right diagram in (j2.5p . A string 

(a, l'f\f'^i^^ is singular if = namely if the rigging attains the allowed maximum ()2.17p . 
We also recall that one actually has to attach the data L with (^, r) to specify the vacancy 
numbers pf^ ()2.7p . Thus we write a rigged configuration as (/x,r)i. We regard a highest path 
6i <8) • • • ® 6l, G (iji.i)®^ as a word 6162 • • • £ {li 2, . . . , n + 1}^. The algorithms explained 
below obviously satisfy the property (j2.20p . 
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C.2 Algorithm for (f> 

Given a highest path bi . . . bi, we construct the rigged configuration (j){bi . . . 6l) = (/x, r)^, in- 
ductively with respect to L. When L = 0, we understand that 0(-) is the array of empty 
partitions. Suppose that 0(6i . . . bi) = {fj,, r)i has been obtained. Denote G {1, . . . , n + 1} 
simply by d. We are to construct {fi',r')i+i = (j){bi . . .bid) from {fJ-,r)i and d. If d = 1, then 
{fi',r')L+i = (/i,r)i+i meaning that nothing needs to be done except increasing the vacancy 
numbers pf'^ by one. (Recall in ()2.7p that only p^^^ depends on L.) Suppose d > 2. 

(i) Set i^'^^ = oo. For c = d — l,d — 2, . . . ,1 in this order, proceed as follows. Find the color c 
singular string whose length i^^^ is largest within the condition i^'^^ < If there are 

more than one such strings, pick any one of them. If there is no such string with color c, 
set i^'^'^ = 0. Denote these selected strings by (^c,£^'^\ri^^^ with c = d — l,d — 2, . . . , 1, 
where it is actually void when i^"^^ = 0. 

{ii) Replace the selected string (^c,£^'^\ri^^) by (c, ^^^^^ + l,ri'^^) for all c = (i — 1, d — 2, . . . , 1 
leaving the other strings unchanged. Here the new rigging r, is to be fixed so that 
the enlarged string (^c,£^'^^ + becomes singular with respect to the resulting new 

configuration fi' and L + 1. 

The algorithm is known to be well-defined and the resulting object gives the sought rigged 
configuration (/i',r')L+i = (p{bi . . . bid). 



C.3 Algorithm for </> ^ 

Given a rigged configuration (/i, r)i , we construct a highest path bi . . .bi = ((/-f, r) l) induc- 
tively with respect to L. We are to determine d{= bi) G {1, . . . , n + 1} and (/i', r') i_i such that 

4>-\{^l,v)L) = <t^-\{^l',v')L-l)d. 

ii) Set = 1. For c = 1,2, ... , n in this order, proceed as follows until stopped. Find the 
color c singular string whose length l^^'^ is smallest within the condition < If 

there are more than one such strings, pick any one of them. If there is no such string with 
color c, set d = c and stop. If c = n and such a color n string still exists, set d = n + 1 
and stop. Denote these selected strings by (c, £^'^^ , ri^^) with c = 1,2, . . . , d — I. 

(ii) Replace the selected string (c, I'^'^^ri^'^) by (c, £^'^^ — l,ri^^) for all c = 1, 2, . . . , d— 1 leaving 
the other strings unchanged. When £^^^ = 1, this means that the length one string is to be 
eliminated. The new rigging ri'^^ is to be fixed so that the shortened string (c, £^'^'^ — 1, ri'^^) 
becomes singular with respect to the resulting new configuration fi' and L — 1. 

For empty rigged configuration, we understand that (j)^^{{0, 0)l) = '/'"^((^i 1 = • • • = 

L 

11 ... 1. The algorithm is known well-defined and ends up with the empty rigged configuration 
at L = 0. The resulting sequence gives the sought highest path 6i . . . 6^ = (/)~^((^, r)/,). 

Example C.l. Let us demonstrate along an L = 8 example. 



1 






1 1 


CO 




2 








1 





11213122. 



(C.l) 
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This is the bottom right case of Example 12.41 The numbers on the left of Young diagrams are 
vacancy numbers, which are subsidiary data but convenient to keep track of. 



L = 8 



L = 7 



L = 6 



1 






1 1 


3 




2 








1 





= 5 

oD 



oD 



L = 4 
2^1 



L = 3 

1 Qi 



L 









Here we have exhibited the number d in the algorithm (i) over the arrows. After this, we only get 

22131211 

1 twice. Thus we get the sequence i— t-i— t-i— t-i— t-i— t-i— t-i— t-i— T-to reach the trivial rigged configuration. 
Reading these numbers backwards we obtain (jC.ip . 

D Proof of Theorem [23] 



For p = bi 
diagram 



G B, define another version of energy E^^\p) 



by taking v 



the highest element of 



whose entries in the jth row are all j. For example, u 



2,3 



111 

222' 



ei + ■ ■ ■ + cl with the 
u^'' is the semistandard tableau 
Set 



V{fi) = {p^V\E\''\p) 



aa 



min ( Z , / • ) m • 



i=l 



(D.l) 
(D.2) 



■p+(/x) = {p G 'P(/i) I p : highest}. 

These definitions resemble (j2.1ip and ()2.12|) . The differences are that here enters E^l^^ instead 
of E^^\ and there is no requirement that p be evolvable in (jP.ip . 

Theorem D.l ([33])' The KKR bijection (j2.18p can be restricted into a finer bijection: 
Lemma D.l. The equality e\^\p) = e\^\p) is valid for any p E Vj^{^j). 

Proof. Take p G V+{^). Then by Lemma 14.11 the rigged configuration for p ® p has the 
configuration obtained from /i by replacing m|"^ with 2nrif^ . Thus from (jP.ip and Theorem lD.il 



one has E^ [p®p) = 2Ej' [p). On the other hand, define u'"'' £ B^'^ (and p) by u^'^^p ~ p^v"^ 
Then, using Lemma 14.11 in this paper and Lemma 8.5 in [18], one can show u'"'' (E> p — p' (E> v' 
for some p' £ B. Since p is evolvable, we can claim that f^'' works as a carrier to define E]j^\p). 
Then from the relation u^'' i^i p (d) p — p 'S' f*"'' p — p 'S' p' 'S' f^'', the energy E^^^ associated to 
p ® p is given by Ei''\p (i^ p) = e\^^ {p) + E^''\p). Comparing this with E^''\p p) = 2Ej^^\p), 
we get e\'^\p) = e\''\p). ■ 

Lemma D.2. There is a natural inclusion V-^{fi). 



Proof. Due to Lemma |D.1|, the difference of V+{fJ.) (|D.2p and V+{fi) (|2.12p is only the extra 
condition that p be evolvable in (j2.1ip . ■ 
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We expect that actually V+ifJ-) = "P+ifJ-) holds, namely, all the highest paths are evolvable. 
However we do not need this fact in this paper. 

Proof of Theorem 12.21 Combine Theorem ID . 1 1 and Lemma [0.21 ■ 
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